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1. INTRODUCTION 


Iris the aim of this address to give a fresh orientation to the subject of crystal 
physics. We begin with a re-examination of the fundamental notions on 
which a theoretical approach to the problems of the subject has to be based. 
The ideas which emerge are then applied to find an answer to the following 
questions: What is the nature of the thermal agitation in crystals ? What is 
the character of their vibration spectra ? These questions are of fundamental 
importance since they confront us in any attempt to explain or interpret the 
physical behaviour of crystals. The results which emerge from our enquiry 
impinge on many aspects of crystal physics. But we shall restrict ourselves 
here to a consideration of those in which the relationships between theoretical 
reasoning and factual experience are most evident, viz., the thermal properties 
of crystals and their spectroscopic behaviour. Even in these fields, we can 
here do no more than broadly indicate the general nature of the results which 
follow from the application of the new ideas. 


2. THE THEORY OF EINSTEIN 


It is appropriate that we begin with a reference to the classic paper of 
1907 by Einstein introducing the quantum theory of specific heats, a paper 
which even fifty years after its publication remains worthy of careful study. 
Einstein commences his paper by showing that the acceptance of Planck’s 
fadiation formula demands as a necessary consequence that the energy of 
a resonator capable of absorbing or emitting radiation with a specific fre- 
quency can only increase or diminish by quanta proportional to that fre- 
quency. It is inferred that the energy of any elementary structure capable of 
mechanical vibration should likewise obey the quantum rule. Einstein then 
indicates a procedure which enables this principle to be applied to the case 
of crystals. He identifies the thermal energy of a crystal with the mechanical 
energy of vibration of an immense number of individual oscillators distri- 
buted over its volume, the total number of these oscillators being the same 
as thrice the number of atoms comprised in the crystal. It is assumed that 
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the oscillators could be grouped into sets; each set comprises a great number 
of individual oscillators characterised by a common frequency of vibration. 
The individual oscillators in each set occupy the various energy levels allowed 
by the quantum hypothesis for the particular frequency. The numbers 
occupying these levels are in the relative proportions determined by the 
respective energies and by the temperature as indicated by the principle of 
Boltzmann. The average energy of an oscillator in each of the sets is then 
evaluated. Summing up over all the oscillators in each of the sets, and then 
over all the sets, Einstein obtains an explicit formula for the thermal energy 
of the crystal as a function of the temperature. 


Examining the procedure thus summarised, it is evident that every step 
in the argument can be justified. A crystal is an assembly of a great num- 
ber of similar and similarly situated groups of atoms, the normal modes 
and frequencies of vibration of which would necessarily all be identical. 
The total number of degrees of dynamical freedom of the entire system may 
accordingly be divided up into sets, each set representing a great number 
of individual oscillators having a common frequency. Since these oscillators 
are interconnected, they can exchange energy with each other and thereby 
constitute a system in thermodynamic equilibrium of which the behaviour 
could be statistically described in terms of Boltzmann’s principle. In other 
words, the procedure followed effects a synthesis of the results of classical 
dynamics with the notions of the quantum theory and the basic principles of 
thermodynamics. The simplicity of Einstein’s theory is as noteworthy as 
its success in explaining the general character of the specific heat-temperature 
curve of crystals. 


We conclude, therefore, that Einstein’s approach to the problem is 
fundamentally correct. His theory indicates that the atomic vibration spec- 
trum of a crystal would consist of a set of discrete monochromatic frequencies, 
their number being the same as the number of sets of oscillators which are the 
carriers of the thermal energy. By virtue of Boltzmann’s principle, the 
distribution of this energy over the volume of the crystal would exhibit 
fluctuations both in space and in time, the magnitude of the fluctuations 
being the greater the higher the frequency of the oscillators under considera- 
tion. However, the theory as it stands gives no precise indication of how 
the oscillators in each set are to be enumerated and how their respective 
frequencies of vibration are to be determined. Neither is any indication 
given whether the modes of vibration of the solid recognised by the classical 
theory of elasticity have any place in the picture. We shall now proceed 
to show how these lacune in Einstein’s theory may be filled up. 
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3. THE DYNAMICS OF CRYSTAL LATTICES 


The fundamental role which the enumeration of the degrees of freedom 
of atomic movement plays in specific-heat theory derives its justification 
from the well-known theorem in classical mechanics regarding the small 
vibrations of a connected system of particles about their positions of 
equilibrium. This theorem states that all the possible oscillations of the 
system are superpositions of its normal modes of vibration in each of which 
all the particles of the system vibrate with the same frequency and in the 
same or opposite phases. On the other hand, in Einstein’s theory we are 
concerned with the modes and frequencies of vibration of the individual 
oscillators of which the total number is equal to the number of degrees of 
dynamical freedom of the entire system. These two points of view can be 
reconciled by identifying the normal modes and frequencies of vibration 
which are characteristic of the structure of the crystal with those that describe 
the dynamical behaviour of the elementary units of which it is composed. 


The fundamental property of crystal structure is that it comes into 
coincidence with itself following a unit translation along any one of the 
three axes of the lattice. Hence, the normal modes of vibration charac- 
teristic of the structure of a crystal should satisfy a similar requirement. 
This can evidently happen in two ways, viz., the amplitudes and phases of 
oscillation of equivalent atoms in adjacent cells are the same; alternatively, 
the amplitudes are the same while the phases are all reversed following the 
unit translation. Since these two possibilities exist for each of the three 
axes of the lattice, we have 2x22 or eight possible situations. In each 
of these situations, the equations of motion of the p atoms contained in the 
unit cell can be completely solved, yielding us 3p solutions. Thus in all, 
we have 24p solutions from which the three simple translations must be 
excluded. We are then left with (24 p-3) normal modes and frequencies 
of vibration. Thus, the fundamental result emerges that a crystal consist- 
ing of p interpenetrating Bravais lattices of atoms has (24 p—3) characteristic 
modes of vibration, each of which is characterised by a specific frequency. 
In (3 p-3) of these modes, equivalent atoms have the same amplitudes and 
phases of oscillation in the adjacent cells, while in the 21 p other modes the 
amplitudes are the same, while the phases alternate in adjacent cells along 
one, two or all three of the axes of the lattice. 


The (24 p-3) normal modes of vibration indicated by the precedin; 
argument may obviously be regarded as the modes of internal vibration of 
the group of 8 p atoms comprised in a supercell of the crystal lattice whose 
linear dimensions are twice as large as that of the unit cell containing P 
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atoms. The 3 omitted degrees of freedom would then represent the 3 degrees 
of translatory freedom of movement of the whole group of 8 p atoms included 
in the supercell. Thus, it emerges that the structural unit whose dynamical 
behaviour isfrepresentative of the entire crystal is not the unit cell of the 
crystal structure but is twice as large in each direction as the latter. 


4. THE VIBRATION SPECTRUM OF A CRYSTAL 


The considerations set forth above indicate that a crystal may be re- 
garded as an assembly consisting of (24 p-3) sets of similar oscillators, each 
with its own distinctive vibrational frequency. Even a monatomic crystal 
for which p = 1 would have 21 characteristic frequencies, while a diatomic 
crystal for which p = 2 would have 45 and so forth. The number of dis- 
tinct frequencies would however be reduced considerably if the crystal belongs 
to a class with a high degree of symmetry. For instance, it can be readily 
shown that a face-centred cubic lattice of atoms would have only four charac- 
teristic frequencies of vibration, while cubic crystals such as diamond and 
rock-salt having two non-equivalent atoms in the unit cell would each have 
nine distinctive frequencies. The important point to be emphasised is that 
the vibration spectrum of a crystal consists of a set of discrete monochromatic 
frequencies whose number is the greater, the more numerous the atoms 
included in its unit cell are or the lower the symmetry of the crystal. 


To complete the picture, we have to consider the nature of the move- 
ments associated with the three omitted degrees of freedom representing the 
translatory movements of the supercell of the crystal lattice. Such transla- 
tions would necessarily result in local stresses and strains of the same 
general nature as those contemplated in the theory of elasticity and hence 
would also give rise to localised oscillatory movements. In other words, 
to complete the picture of the thermal agitation envisaged in Einstein’s 
theory, we have to superpose on the oscillations of precisely defined fre- 
quencies determined by the atomic structure of the crystal, other quantized 
oscillations whose frequencies are not so precisely defined, since they would 
depend upon the dimensions of the moving masses and the forces thereby 
brought into play. It is evident that the smaller the mass of the oscillators 
under consideration, the more numerous they would be and also the higher 
their frequencies of vibration. It follows that the spectrum of the thermal 
agitation represented by such movements would be continuous and would 
exhibit a concentration of frequencies most marked at the upper end of the 
frequency range and falling off rapidly to zero as we move towards the lower 
end of the range, 
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5. THE THEORIES OF DEBYE AND BORN 


The views set forth above regarding the nature of the thermal agitation 
in crystals and the character of their atomic vibration spectra are essentially 
a justification and a development of the original ideas of Einstein outlined 
in his fundamental paper of 1907. It is a matter of scientific history, how- 
ever, that an approach of an entirely different nature to the problems of 
crystal physics was put forward a few years later and found general acceptance. 
The reference here is to Debye’s theory of specific heats familiar to all stu- 
dents of physics and chemistry and to the lattice dynamics of Born and his 
collaborators which is regarded with favour by the mathematical theorists. 
The theories of Debye and Born differ in detail, but the basic idea in both 
is the same, namely that we are here dealing with a boundary value problem 
in the theory of wave-propagation. Actually, what we are concerned with 
is the determination of the modes and frequencies of vibration of the atoms 
in a crystal about their positions of equilibrium. These modes and fre- 
quencies would be determined by the masses of the atoms and their inter- 
actions with the surrounding atoms in a domain the extent of which is limited 
by the range of the interatomic forces. It follows that considerations re- 
garding wave-propagation and boundary conditions are not relevant to the 
problem. 


A rigorous mathematical investigation of the theory of the propagation 
of waves in a general crystal consisting of p interpenetrating lattices of atoms 
on the basis of the classical mechanics shows that only in the limiting cases 
of the lowest frequencies of vibration is the group velocity of the waves 
identifiable with their phase velocity. With increasing frequency, the phase 
velocity and the group velocity diverge from each other, and the group 
velocity becomes zero for (24p-3) distinct frequencies, the nature of the 
atomic modes of vibration for these frequencies being identical to the 
(24 p-3) normal modes of vibration of the structure of the crystal referred to 
earlier in this address. Thus, the investigation from the standpoint of the 
theory of wave-propagation and of classical mechanics leads to the same result, 
namely that these are the (24 p-3) characteristic modes and frequencies of 
vibration of the structural units of which the crystal is composed. 


We may put the matter a little differently by saying that while Einstein’s 
theory identifies the thermal energy of a crystal with the quantized vibrational 
energy of immense numbers of oscillators having identical frequencies, Born’s 
lattice dynamics identifies it with the quantized energy of immense numbers 
of wave-motions whose frequencies are all different. The two approaches 
are clearly irreconcilable, since in the one case, the vibration spectrum of 
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the crystal consists of a set of discrete frequencies, while in the other, it is 
a continuous spectrum. That Born’s approach to the specific heat problem 
is indefensible becomes evident when we remark that his “ waves” are not 
normal modes and hence are incapable of enumeration. 


6. THE THERMAL PROPERTIES OF CRYSTALS 


We shall now consider the application of the ideas indicated in the 
present address to the interpretation of the physical behaviour of crystals. 
The thermal properties may be considered first, since the nature of the thermal 
agitation plays a leading role in relation to them. 


A familiar fact of experience is the extreme slowness with which heat 
diffuses through a solid. On the other hand, elastic waves travel through a 
solid with a velocity of some thousands of metres per second. It would 
be a formidable problem to explain such a fundamental difference in beha- 
viour, if one were to identify thermal energy with wave-motion, as is done in 
the theories of Debye and Born. Such a paradox does not however con- 
front us on the present view of the nature of the thermal agitation in which 
it is identified with the energy of vibration of individual oscillators localised 
in the crystal. The conduction of heat through a solid would on this view 
be the result of an exchange of energy quanta between the individual oscilla- 
tors. Such a process would bear no resemblance whatever to wave-propa- 
gation in the ordinary sense of the word. 


The evaluation of the specific heats of crystals as a function of tempe- 
rature reduces itself to the problem of determining the (24 p—3) characteristic 
frequencies of atomic vibration. When p is unity or a small integer, the 21 p 
superlattice frequencies make important contributions to the thermal energy. 
When p is large, however, the (3 p-3) modes and frequencies of the unit cell 
reappear with only minor changes as superlattice modes and frequencies, 
and hence many of the latter could be identified with the former without 
sensible error in the calculation of the thermal energy; a considerable simpli- 
fication thereby results. The thermal energy represented by the three trans- 
lations of the supercell may be assumed to be proportional to the absolute 
temperature, except at the very lowest temperatures where its continuous 
spectrum of frequencies would have to be explicitly to be taken into consi- 
deration. 


Crystals exhibit a surprising variety of behaviour in respect of the 
changes in their linear dimensions with rise or fall of temperature. The 
present view of the nature of thermal agitation which regards it as a localised 
vibration of individual oscillators opens a pathway to the interpretation of 
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the observed facts. The finite amplitudes of vibration of the individual 
oscillators and the mechanical anharmonicity associated therewith would 
result in the frequencies of vibration diminishing in the higher states of 
excitation. It would also result in a change of the average dimensions of 
the individual oscillators, which again would depend on the modes of 
vibration and hence may be very different in different directions. The 
integrated effect would be a shifting accompanied by a broadening of the 
lines in the vibration spectra and a change in the linear dimensions of the 
crystal. Thus, a close correlation is to be expected between the thermal 
expansion of crystals and the changes in their spectroscopic behaviour with 
rise or fall of temperature. 


7. THE SPECTROSCOPIC BEHAVIOUR OF CRYSTALS 


The vibration spectra of crystals are accessible to experimental study 
with a great variety of materials and by diverse techniques of investigation. 
The nature of the spectra becomes manifest when the absorption of light 
by crystals exhibiting visible colour (e.g., the inorganic chromates), or of 
ultra-violet light by crystals (e.g., the aromatic organic compounds) is in- 
vestigated with the crystals cooled down to the lowest possible temperatures. 
The absorption spectra then appear fully resolved into a set of sharply 
' defined lines. These can be interpreted as combinations of the electronic 
with the vibrational frequencies of the crystal. Very low temperatures are 
also necessary for effecting a similar resolution into sharp lines of the lumi- 
nescence spectra of various crystals of which such luminescence is an inherent 
property, ¢.g., the uranyl salts and of other crystals of which the luminescence 
is ascribable to the presence of an impurity phosphor (e.g., ruby). Such 
luminescent crystals also exhibit at low temperatures characteristic absorp- 
tion spectra which are resolved into discrete lines and which are related by 
an approximate mirror image symmetry to their emission spectra. 


There are also other and rather more direct methods of investigation. 
Three different techniques based on infra-red spectroscopy are available, 
viz., the emission, absorption and reflection of infra-red radiations by crystals. 
In these techniques, the disturbing effects of thermal agitation and of the 
anharmonicity of the vibrating electric doublets in the crystal are strongly 
manifested. Nevertheless, by an appropriate choice of the material to be 
investigated and of the techniques used, e.g., by varying the thickness of 
the crystal, cooling it to low temperatures and the application of adequate 
resolving power, it is possible to obtain results capable of unequivocal inter- 
pretation. Finally, we have the method based on the study of the scattering 
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of light in crystals. The results obtained by this method are of such im- 
portance that they will be considered separately in the next section. 


The whole body of experimental evidence forthcoming with the mate- 
rials and methods described above confirms the thesis that apart from dis- 
turbing effects such as those mentioned earlier, the atomic vibrations in 
crystals appear as a set of sharply defined monochromatic frequencies in 
the spectrum. 


8. THE SCATTERING OF LIGHT IN CRYSTALS 


The nature of the results to be expected when a crystal is traversed by 
an intense monochromatic beam of light is different in relation to the (3 p-3) 
modes of atomic vibration which have the same phase in adjoining cells of 
the lattice, and the 21 p modes in which the phases are opposite in adjoining 
cells along one, two or all three of the axes of the lattice. In either case, 
we are concerned with the scattered radiations having their origin in the 
changes of optical polarisability of the vibrating atomic groups. Ordi- 
narily, except in the cases where the symmetry of the vibrational 
modes results in there being no resultant change of optical polarisability, 
the (3 p-3) modes would manifest themselves as frequency shifts in the 
spectrum of the scattered light. On the other hand in the 21 p modes, the 
optical moments being in opposite phases in adjoining cells, their effects 
would cancel out. However, if the amplitudes of vibration are sufficient, 
the changes in optical polarisability proportional to the squares of the atomic 
displacements would also have to be considered. It is then possible for all 
of the (24 p-3) modes of atomic vibration to appear with doubled frequency 
shifts in the spectrum of the scattered light. The intensity of such scattering, 
being proportional to the square of the absolute atomic displacements, would 
necessarily be small. Strongly exposed spectrograms would accordingly 
be necessary to record the same. 


It is evident that the study of the second-order spectra of light- 
scattering is of the greatest importance, since in the first place they allow 
us to observe and measure the frequencies of modes which do not appear 
in the first-order spectrum, and secondly because the intensity of the lines 
of the second-order spectra relative to those of the first-order spectra is a 
measure of the absolute magnitude of the interatomic displacements res- 
ponsible for the changes in optical polarisability. Studies of this kind are 
obviously most suitably undertaken with crystals of relatively simple struc- 
tures, as for example, diamond, rock-salt, calcite and quartz, since the 
results are then most readily interpreted. 
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Very detailed investigations of the kind indicated carried out at 
Bangalore have established that the vibration spectra of crystals include 
both the (3 p-3) and the 21 p modes with discrete monochromatic fre- 
quencies, and further that the oscillators whose vibrations manifest them- 
selves as frequency shifts are of the same order of magnitude in their dimen- 
sions as the unit cells of the crystal structure. The latter result is inferrable 
from the fact that the relative intensity of the first- and the second-order 
spectra is of the same order of magnitude in crystals as in fluids. For, the 
size of the quantized oscillators emitting the scattered radiation of any 
given frequency determines the absolute magnitude of the interatomic dis- 
placements responsible for the variations of optical polarisability. 


SUMMARY 


1. Einstein’s view of a crystal as an assembly of immense numbers of 
quantized oscillators having a common set of vibration frequencies is not 
only the logical and correct view of the matter but also proves itself when 
fully developed to be an eminently successful view. It gives us a deep and 
quantitative insight into the thermal behaviour of solids. The nature of 
the vibration spectra of crystals indicated by it is confirmed in detail by the 
results of spectroscopic studies with diverse materials and by various 
techniques and with especial completeness by studies on the diffusion of 
light in crystals. 


2. The ideas underlying the specific-heat theories of Debye and Born 
are irreconcilable with the observed spectroscopic behaviour of crystals and 
especially with the effects exhibited in their second-order spectra of light- 
scattering. 
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I. INTRODUCTION 


IN a recent paper,! observations of the intensity and polarisation of the 
light from the sky during twilight taken at Mount Abu during 1951 were 
presented and briefly discussed. It was shown that for depressions of the 
sun exceeding 7° or 8°, the zenith sky was much brighter than would be 
expected if primary scattering alone were responsible for the light and that 
in the treatment of twilight problems, due account should be taken of the 
contributions from secondary scattering, as has been shown by Hulburt.? 
Marked discontinuities in the polarisation of the light from the zenith sky 
at solar depressions of 5°, 8°, 12° and 17° support the presence in twilight 
of other contributions in addition to primary scattered light. 


1.1. Summary of previous work 


Robley* has recently given a method of calculating the intensity and 
polarisation of the secondary (multiply) scattered light during twilight. His 
work is based on the work of Dr. S. Chandrasekhar* for an infinite plane- 
parallel atmosphere of defined optical thickness. To compute the secondary 
scattered radiation for the zenith sky for a solar depression D, Robley extends 
the atmospheric path traversed by the solar rays to a vertical plane in the 
atmosphere passing through the centre of the earth and making an angle 
D/2 with the observer’s vertical. He then calculates the scattered light 
emanating from this plane in a direction parallel to the observer’s horizon, 
with certain simplifying assumptions which lead to a constant value of 50% 
for the polarisation of the secondary scattered light from the zenith sky 
throughout the period of twilight. Prima facie, this does not appear to be 
correct, because observation shows that with increasing values of solar depres- 
sion, the brightness of the sky becomes more and more concentrated towards 
the horizon in the sun’s meridian plane. 


1.2. General outline of the method adopted in this paper 


The method adopted in this paper for the calculation of secondary 
scattering is based on the work of Chapman and Hammad.5:* The following 
assumptions were made :— 
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(1) The earth was taken to ‘be a uniform sphere of radius 6370 km., 
surrounded by an exponential atmosphere of equivalent scale 
height 7 km. 


(2) The atmosphere was arbitrarily divided into a number of concentric 
shells of 4km. thickness, for computing the intensities of prim- 
ary and secondary scattered radiations. 


(3) The solar radiation falling on the top of the atmosphere was assumed 
to be unpolarised and its intensity taken to be unity. 


(4) The depolarisation factor d of air due to molecular anisotropy was 
taken to be 0-04. 


(5) Higher orders of scattering than secondary were neglected. 


(6) No correction was made for the absorption or scattering of light 
by dust or haze. 


(7) Atmospheric refraction was neglected. 
(8) No correction was made for reflection by the ground. 


If. THEORY 


In Fig. 1, ABC represents a section of the earth; O is its centre and a 
its radius. R, is the position of the observer on the surface of the earth. 
OR, represents the direction of the observer’s zenith. Q is a point in the 
atmosphere in the direction of observation R,Q. The Z-axis is chosen to 
lie along the observer’s zenith OR,. The X-axis lies in the plane of the 
paper such that the XOZ plane is parallel to the incident solar radiation. 
The Y-axis is perpendicular to the plane of paper passing through O. The 
polar co-ordinates (8, ¢) of a given direction are such that @ is the angle 
which the direction makes with the Z-axis, and ¢ is the azimuth with XOZ 
as the reference plane. With this convention, the polar co-ordinates of 
the incident solar radiation SP are (Z,, O). If the direction of observation 
is the zenith or R,Q, its 0, ¢ are (O, O). PQ the direction of primary scat- 
tered radiation is given by (6’, ¢’). The unit vectors along SP, QR, and PQ 
are denoted respectively by s, k and k’. 


2.1. Intensity of primary scattered radiation 


Consider a solar beam of intensity I,, travelling along SP. I, denotes 
the intensity of solar radiation outside the earth’s atmosphere in a specified 
narrow wave band, being the energy received per unit time on unit area normal 
to the vector s. Let o denote the mass coefficient of scattering and p the 
density of air at a distance s from P along PS. 














O 


Fic. 1. The incident light is primarily scattered at P (where the local zenith angle of the 
sun is Z, and the height R,P = c,,). The direction PQ is defined by (6’, g’). The direction of 
observation towards the observer’s zenith is defined by @ = 0, ¢ = 0. The primary scattered light 
travelling along PQ gets rescattered at Q along QR, towards the observer. QR, = J, 


Then 
P 
[o-p-ds = Tsp (1) 


is the equivalent optical path which the solar radiation has to travel in the 
earth’s atmosphere before reaching P. It will be shown later that rsp is a 
function of the vertical height cm of P above the surface of the earth and 
of the local zenith angle (Z) of the sun at P. The intensity of the solar radia- 
tion at P is then given by 


I=I.e7’s (2) 


Let PQ (6’, 4’) be a direction along which the incident solar radiation 
is scattered from P according to the Rayleigh law of scattering. The energy 
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scattered from P in a cone of small solid angle dk’ with apex at P is given 
by 
men 3ol,,, 
Pid = ton (+ #2) 
where d is the depolarisation factor due to molecular anisotropy and yf’ the 
angle between s and k’ given by 


cos %s’ = cos Z, cos 6’ + sin Z, sin 9 cos ¢’ (4) 


e~Tse (1+cos*b’+d sin? %’) dk’ = (3) 


This scattered radiation is polarised and can be resolved into two beams 
of plane polarised light at right angles to each other. If Ej», represents 
the primary scattered emission from P in the direction of k’ with its electric 
vibration along the vector n’ perpendicular to the plane sk’ at P, and E,z,, 
the intensity of the radiation with its electric vector along ¢’ perpendicular 
to n', we have 


Eww = Tee + Ha) ©" 6) 
and 

Ew = KE, (6) 
where 

K = d+ (1 — @) cos*’ (7) 


Consider a small surface area dS at Q perpendicular to PQ. The 
amount of primary scattered radiation from P passing through dS per second 
is given by 

Ein (1 + K) e-?re X 2) 
where dS/s? is the solid angle subtended by dS at P and 7p, the equivalent 
optical path from P to Q. It will be later shown that tpg which is given 
by 


Q 
Tp = J o-p-ds 
P 


is a function of 6’, cm and by. If dV isa small volume element surrounding 
P and subtending a small solid angle dk’ at dS, and if p’ is the density at P, 
then the amount of primary scattered radiation passing through dS per 
second from a mass element p’.dV surrounding P is given by 


Ein (1 -+ K) e-tm & p'.dV 
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and the intensities R,,, and R,z of the normal and transverse components 
of primary scattered radiation passing normally through unit area at Q 
through unit solid angle along k’ is given by 


<2 dS, p-dV 
w =f Ean ee X OS dk 


and 
Ry = KRiv 


the integration extending from co to Q when the whole path is illuminated 
by solar rays, or upto the shadow limit when a part of PQ lies in the shadow 
region of the earth. 


Substituting the values of E,,, and of dV = s*.ds.dk’ in Rin, we have 





e+ + Fd) ¢ fe e-7 see— Trop’ -ds (8) 
and 
Rit = KRin (9) 


It must be noted that unlike the plane-parallel atmosphere, cos Z on 
the actual earth changes along QP and. also with ¢’. Hence 7p is not inde- 
pendent of ¢’. Thus the integration formule developed by Chapman and 
Hammad for Ryn, are not applicable in this and we have adopted numerical 
integration. 


2.2. Intensity of secondary scattered radiation 


Having obtained the values of Ryn, and Rj, for different values of 6’ 
and ¢’ at Q, the next problem is to consider the effective contribution by 
the volume element at Q towards the secondary scattered radiation travelling 
along QRy. Let m be the unit vector perpendicular to the k’k plane at Q. 
If Eyn and E,; are defined as the energy of the secondary emission for normal 
and transverse components emitted at Q respectively, then as shown by 
Hammad, we have 


Ean = 16m (1 Be 3d) / {Ay 2 KA3} Rin dk 


and 


30 
Tor (i + $d) J (Br + KBs} Ran dk 


Ext = 
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where 


A, = d+ 2(1 — d) cos? Xnz 
B, = d+ 2(1 — d) cos? Xtp, (10) 
B, = d + 2(1 — d) cos? Xe, 


A, = d+ 2(1 — d) cos? Xnn, 


Xnn> Xnt» Xtnr and X¢tt, are respectively the angles between the vectors n and 
n',nand t', t and n’ and ¢ and ¢’ and the formule for obtaining their values 
are as given by Hammad. 


COS Xnn = cos (€ — €’) (11) 
COS Xnt = cos #’ sin (€ — £’) (12) 
COS Xtn, = cos % sin (€’ — &) (13) 


cos Xz = sin ¢% sin %’ + cos % cos #’ cos (é' — &) (14) 
where #, v’, € and é’ are defined by the equations 


cos %* = cos Z, cos 8 + sin Zp, sin 8 cos ¢ (15) 
and 
. . Sin 6 sind 
sin € = inf (16) 


Since dk’ = sin 0’.d6’.d’, the equations for E,, and E, can be re- 
written as 


Een = 7¢- lord yay JJ (A+ KAg} Ran sin 6°.d6"-d9" (17) 


and 
3 hs Ul , , 
Est = lee +4 / Jf {B:+KB,} Ruy sin 6’.d6’.d¢’ (18) 


Knowing Eg, and E.¢ for a number of points along QR,, the tota 
secondary scattered radiation received at R, along QRo, can be integrated. 
Its normal and transverse components Re, and R,z are given by 


Ran = / Bane 7 or p-ds (19) 
Ro 
and 
Rat = / Eot €~70% p.ds (20) 
Ro 


where p is the density of air at Q and 79x, is the equivalent optical path from 
Qto Ry. It is a function of by, as well as of 0. 











342 J. V. DAVE 


Ill. TABLES FOR COMPUTING THE SECONDARY SCATTERED RADIATION 


The quantities required for computing R,», by the summation method 
such as s, 4s, tsp, have been tabulated by many workers such as Sekera,’ 
Chapman®.® and others. However, their tabulation could not easily fit into 
our system, and all the quantities required for this work were therefore 
tabulated from the formule given below. 


As already mentioned in 1.2, the atmosphere is arbitrarily divided into 
a number of concentric layers of 4 km. thickness. The base of the nth layer 
is at a height of 4(m — 1) km. while its centre is at (4m —2)km. The 
positions of Q are selected at the base of the different layers, and of P at 
the centres of the successive layers. On replacing the integral of Equation (8) 
by a sum, we have 


3opol.., sp e7T. c 
Riw = [6 (1 + 4d) P 2 ilies alia (21) 


p’ = pye~°m/* for an exponential atmosphere and H the equivalent scale height 
is assumed to have a uniform value of 7 km. 


If the values of rs), tpg at the centres of the different layers and also 
of 4s the thickness of the layers for any 6’ are known, Ryn, can be calculated. 
Since tcp is a function of Z and cm, we have to calculate for a given position 
of Q and a fixed value of Z,, the local zenith angle Z of the sun at the centres 
of the successive layers along QP. For this, it is necessary to tabulate a 
shown in Fig. 1 for different values of Cm. 


3.1. s, the distance of the centre of the layer from the base of another layer 
in any direction 


From Fig. 1, we have 
(a + Cm)? = s? + (a + by)? + 285 (a + dy) cos 0’ 
On neglecting cm? and 5b,” in comparison with the other terms, 
s? + 2s (a + by) cos 0’ — 2a(cm — bn) = 0 (22) 


With the help of Equation (22), the values of s were tabulated for differ- 
ent values of cm and 6’ = 0°, 15°, 30°, 45°, 60°, 75°, 80°, 85°, 88° and 90° 
when the position of Q is at b,, the ground level. From this table, it is easy 
to obtain the values of s corresponding to any b, and 6’ between 0° and 90° 
with sufficient accuracy. For positions of Q above the ground, when 0’ 
can take values greater than 90°, QP can intercept some of the lower layers 
twice. This was appropriately allowed for in computing the path lengths 
for 6’ between 90° and 180° at similar intervals. It is expected that the 
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tabulated values of s are correct within 1 km. when s is small and within 
2km. when it is large. 
3.2. As, thickness of the layer in any direction 

Equation (22) was used to calculate the distance 4s between the bases 


of two consecutive layers for all the above-mentioned values of 6’. Loge 4s 
was tabulated for all values of cm, bn and 6 required in the calculation. 


3.3. Z, the local zenith distance of the sun for any position of P for a given Z, 

In Fig. 1, the vertical through P (PR,O) makes an angle a with the Z- 
axis. It can be represented by (a, ¢’) and Z is given by 

cos Z = cos Z, cos a + sin Z, sin a cos ¢’ (23) 

a can be determined from the relation 

s 
at Cm 

The values of sin a and cos a were tabulated for all the positions men- 
tioned in 3.1 and used for the other similar positions as necessary. 


When Z is greater than 90°, there is a limit below which P would not be 
illuminated by direct sunlight. Limiting values of cos Z for different values 
of Cm can be obtained from 


sin a = 





sin 6’ (24) 


: a 
sin Z,; = eo (25) 
With the help of Equation (24) tables of ‘cos Z were prepared for Z, 
corresponding to 90°, 94°, 98°, 102° and 106° for all the values of 6’ and by, 
and at intervals of 15° for ¢’. 
3.4. tsp, the equivalent optical path from outside the atmosphere to the point 
P which is at a height Cm above ground and for different values of cos Z 
Consider a point P situated at a height cm in its vertical above the 
surface of the earth. Let the incident solar radiation make an angle Z with 
OP (Fig. 2). Let P’ be any point in SP, at a distance s from P and at a 
vertical height 4 from the ground. The density at P’ is given by 


p= py 


and the equivalent optical path 7sp (cos Z, Cm) as defined by (1) is given by 
P 
Tsp (cos Z, Cm) = oPy i e Wit ds (26) 


8a? — (ng? — 1)? | 6 +.) 
a ha oe (27) 
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where 
A = wave-length of the incident radiation 
n, = refractive index of the air 
N, = number of molecules per c.c. at the ground 
Py = density of the air in gm. per c.c. at the ground. 
From Fig. 2, we have, 
(a + h)? = s? + (@ + Cm)? + 25 (@ + Cm) cos Z 
and neglecting c,? and h? compared to the other terms, we have, 

















2 
h= = 4 EP Oe 5 + 
180-Z 
we 
Ca 
Pp 
Cm 
a 
a 
re) ° 
Fic. 2 Fic. 3 


On substituting this in Equation (26) and carrying out the integration, we 
have 
—ry (a +€m)* cos? Z 
Tsp (cos Zz, Cm) = Opp af re twlxe 2aH 


i (a + Cm) cos Z | 
[1- Es ee |) 

The function Erf(y) is tabulated in books of Mathematical Tables at 
regular intervals between 0 and 3 while for y > 3, we can use the standard 
expansion formula. It can be shown that for cos Z = 1 and cm = 0, Equa- 
tion (28) reduces to 
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Tsp (1, 0) = oppH = 7y 
= the vertical optical path (29) 
After substituting op, by 7,/H in Equation (28), the equation was used 
to tabulate t<p (cos Z, O)/7, for a number of intervals of cos Z between 0 
and 1. From this, we can calculate 7,» for any wave-length by selecting 7}. 
The intervals are sufficiently numerous to obtain t¢p/7, for any cos Z within 


the accuracy required. For positive values of cos Z and Cm, e~‘m/* is the 
most important factor and 


Tsp (Cos Z, Cm)/7, = Tsp (Cos Z, O) eCm/F/7, (30) 


For negative values of cos Z, the following method is used. From 
Fig. 3, 

a+h 

a+Cm 

from which hy can be calculated for a given Z and cm. Then 


Tsp (Cos Z, Cm) _ 27sp (O, O) nym 
T1 Ty 


= sin Z 








_ Tsp {cos (180 — Z), Cm} (31) 
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Using Equation (31), tsp (cos Z, Cm)/t, were tabulated at regular intervals 
from zero to the shadow limit for different values of c», from c, to cys. For 
higher values of cm the calculations were carried out whenever necessary, 


3.5. tpg, equivalent optical path between any two points in the atmosphere 
From Figs. 4A and B, we can write 


Tp (cos 9’, Cm to bn) _ Tsp (cos 9’, bn) 
71 mF Ty 








__ Tsp {cos (8 — a), Cm} 
£08 UA) Oa 


(32) 


For negative values of 6’, for which PQ intersects the surface of the earth, 
the calculation of tpg is similar to that indicated in Fig. 4C. 


TV. INTENSITY AND POLARISATION OF THE ZENITH SKY 


In this section the results of the computation of the intensity and polar- 
isation of the primary and secondary radiations from the zenith sky during 
twilight for 7, = 0-2 corresponding to 4 = 4350 A are given. The effect 
of adding a fixed intensity of night air-glow radiation with a polarisation of 
4% as obtained from twilight observations is then discussed. 


4.1. Rin, primary scattered radiation received from any direction for different 
values of Z, 


Rin was calculated with the help of Equation (21). The summation 
extends over all the layers which lie in the sunlit part of the atmosphere. 
I, was taken to be equal to unity, and op, = 7,/H from Equation (29). 
Equation (21) shows that the part e~7re e~m/# As is independent of ¢’ and 
Z,. An auxiliary table of {— tpg — Cm/H + Loge 4s}. for different values 
of b, and 7, = 0-2 was prepared to reduce the work of computation. 


When R,n, was computed for different values of 6’ and ¢’, for Z, = 90° 
and by, = bd, or bs, it was noticed that for a fixed 6’, the graph of Rn, against 
cos ¢’ was nearly a straight line for values of 8’ between 0° to 75° and between 
100° to 180° (Fig. 5). It was slightly curved when 6 lay between 75° and 
95°. Therefore, Rin, was obtained by linear interpolation after computing 
its value for cos ¢’ = 1-0, 0-5, 0-0, — 0-5 and — 1-0 for a given by. This 
procedure was adopted when Z, was 90° and 94° and the whole sky more 
or less illuminated by direct solar radiation. For higher values of Zp), Rin 
was calculated individually for each ¢’ and 6’ as only a small portion of the 
sky could then get direct sunshine. 
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Fic. 5. Variation of the normal component of primary scattered light (R,,’) with solar azimuth 
(cos ¢’) for different zenith angles of the sky from 15° to 90° as seen by an observer at ground 
level (6, = Okm.). The sun is assumed to be on the observer’s horizon (Z) = 90°). 

The variation of the normal component of primary scattered radiation 
as received by the observer at the ground from different parts of the sky 


for Z, from 90° to 102° is shown in Fig. 6. 


4.2. Exn and E,4: Secondary emission from any point in the direction of 
observation 


Equation (17) can be rewritten as 


3 ' 
Ean = ry Dy Dy KAD Ra 
sin’. 
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Knowing the values of {A, + KA,} and Rjy, at definite intervals of 6’ and 
¢' for different values of b, and Z,, further calculation was carried out in 
the following manner: 


First, summation was carried out over ¢’, keeping 6’ constant. Ryn’ 
was multiplied by the corresponding {A, + KAg}, the arithmetic means of 
the consecutive products were then added and the sum multiplied by d¢’ 
= (5237 radians corresponding to 30° interval. Symmetry about the 
principal plane gives the value for 30° from the tabulated values for 15° 
intervals. 


For summation over 0’, X[{A, + KA,} Rind¢’] was multiplied by 
the corresponding value of sin 6’, and the arithmetic means of consecutive 
terms were taken. Rjpn, includes the exponential exp(— tsp — Tpg). These 
were then multiplied by corresponding intervals of d6’ in radians. The sum- 
mation of the product multiplied by 30 po/16m7 (1 + 4d) gave the value of 
po Ean for the particular level b, and Z, under consideration. The values 
of po Een and poE.¢ so obtained are given in Table I for different values of 
by, and Zp. 


TABLE I 
py Eon, for the observer's zenith for different values of by and Zy 








Zo 90° - 94° 98° 102° 106° 
b, 

b, 0:245x10-* 0-143x10- 0:015x10-* 0:013x10-* 0-000x10-” 
b, 0-339 0-199 0-032 0-051 0-025 

b, 0-445 0-262 0-062 0-133 0-117 

by 0-535 0-326 0-101 0-287 0-383 

bs 0-584 0-406 0-162 0-605 1-037 

be 0-587 0-491 0-237 1-040 1-994 

by 0-295 1-389 2-769 

bg 0-311 1-507 3-048 


Po E,, for the observer’s zenith for different values of b, and Z, 


by, 0:077x 10 0-044x10-> 0-004x10-* 0-:002x10-® 0-000x10-” 


by 0-113 0-062 0-008 0-008 0-003 
bg 0-156 0-082 0-014 0-021 0-013 
by 0-193 0-105 0-022 0-046 0-045 
bs 0-214 0-132 0-036 0-096 0-125 
be 0-216 0-162 0-053 0-167 0-244 
by 0-066 0-225 0-340 


by 0-071 0-247 0-383 
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VERTICAL DISTRIBUTION PRIM. & SEC. SCATTERED RADIATION 


Fic. 7. Vertical distribution of primary and secondary scattered light from the zenith sky 


during twilight for solar zenith angles from 90° to 106°. 
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It will be seen from Table I that the computation has been stopped at 
b, for 90° and 94°. For heights greater than 5,, the value corresponding 
to b; has been used. A similar remark applies to heights greater than b, 
for 98°, 102° and 106°. This is justifiable because at these heights, the density 
becomes relatively small and the effective contribution to the light received 
by the observer becomes negligible (see Fig. 7). 


4.3. Ren and Ry: Normal and transverse components of the total secondary 
scattered light as received by the observer at the ground from his zenith 


Equation (19) can be rewritten as 
Ron = = poEen e7 7% e-Fm/# As (33) 


the summation extending over all the layers above the point of observation. 
4s = 4km. for 6= 0°, 6 = 0°. 7 ox, is the same as 7pg for 6’ = 0° and 
by, = 6, when the observer is situated at ground level and the direction of 
observation is along his zenith. The value of p)E,, is the arithmetic mean 
of two successive values of b, given in Table I. The vertical distribution 
of primary and secondary scattered radiation from the zenith sky for differ- 
ent solar angles are plotted in Fig. 7. 


The values of Ro» and R,z are given in Table II along with Ry», Ri, 
P, and P, where P,, P, are the percentage polarisations of primary and 
secondary scattered light. P is defined by 


P= 7% , * 100 (34) 


TABLE II 


Intensity and polarisation of the primary and secondary scattered 
radiations from the zenith sky 





Zo Ri, Ry P, % Ro, Ry, P, yA 





90 1°30x10-? 0-052x10-* 92:3 2°-43x10* 0-84x10- 49 
94 6:93x10% 0-31 x10® 91:4 1°52x10° 0-48x10~% 52 
9% 9-0 x10 0-5 x10 89 4-13x10-7 —- 0-94x 10°” 63 
102 a if é 1-26x10-§ 0-20x 10-8 72 
1-95 x 10-° 


0-24 x 10-78 78 
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Fic. 8. Variation of the normal and transverse components of primary and secondary 
scattered light from the zenith sky during twilight. 
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44. Variation of the intensity and polarisation of the zenith sky during twilight 


Figure 8 gives the variation of Ryn, Rut, Ren and Rez with Z). Smooth 
curves were drawn from the values given in Table II and from them, the 
values were read out at each 1° interval to calculate the intensities and polar- 
jsations of the zenith sky at different zenith angles of the sun. 


Figure 9 shows the variation with Z, of (a) the primary scattered 
radiation R,, (b) the primary and secondary scattered radiations together, 
R, + Re, (c) Ri + Re + Ro when Ry is the background night air-glow 
radiation assumed to be constant at 1-92 10-* and (d) the observed intensity 
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Fic. 9. Variation of the intensity of the zenith sky during twilight. 
R, .. Primary scattering alone. 
R, + R, .. Primary plus secondary scattering. 


R, + Rg + Ro.. Ry +R, plus background night air-glow of constant 
intensity. 





of the zenith sky for A = 4,500 A during evening twilight as observed at 
Mount Abu. The magnitude of R, was fixed from observations of the zenith 
sky when the sun was more than 18° below the horizon. 
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It can be seen from the figure that the secondary scattered radiation 
becomes more and more important after 95°. Between 98° and 102°, the 
light received by the observer is mainly secondary scattered radiation. After. 
wards, the night air-glow becomes comparable in intensity to secondary 
scattered light, and by 106°, practically only night sky emission remains as 
the factor contributing to the zenith sky intensity. 


In Fig. 10 are given the curves of polarisation of the zenith sky when 
the different contributions are considered singly or jointly. The curve giving 
the polarisation of the primary scattered light shows that it is 92% at sunset 
and decreases to 88% at 100°. On the contrary, S representing the polar- 
isation of the secondary scattered light shows an increase from 48% at 
sunset to 79% at 106°. Adding the normal and transverse components of 
the two radiations from the zenith sky and calculating the effective polarisa- 
tion, the polarisation gets reduced to 84% at sunset. After 94°, the secon- 
dary scattered radiation becomes relatively more important and the polar- 
isation shows a rapid fall; after 98°, P becomes unimportant and P+§ 
follows S. 

If night sky radiation of intensities Ron = 1-00<10~ and Ro; = 0-92 
x 10-® be added to Ryn + Ren and Ryz + Rez respectively, and the polarisa- 
tion recalculated (Curve P + S + N), the curve gradually changes its course 
when the zenith distance of the sun increases beyond 99°, and the change 
becomes rapid after 102°. When Z is greater than 108°, night conditions 
get fully established. 


When the calculated polarisation curve is compared with that obtained 
by observation, it is seen that while there is agreement in general features, 
there are some important differences. The observed values of polarisation 
are markedly smaller. The transitions are sharper in the observed curve 
than in the calculated one.* 


This is probably due to the following causes :— 


(1) The effect of dust in the atmosphere and of tertiary and higher 
orders of scattering have been neglected. Inclusion of these in 
the calculation will decrease the polarisation. 


(2) The density of the atmosphere has been assumed to decrease expo- 
nentially with height with a scale height of 7km. Actually, there 
are well marked stratifications in the atmosphere and the decrease 
of density with height is not a simple exponential. The difference 


* Recent observations of polarisation at Mount Abu with a better type of polaroid have 
shown 80% polarisation at sunset and 55 to 60% between 8 and 12° solar depression. 
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94° and 98° can be related to the existence of the warm ozone layer 
between 30 and 60km. An attempt is being made to remedy 
these defects in the theory. 


e 
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V. POLARISATION OF THE SKY AT POINTS 30° FROM THE ZENITH AND 
IN THE PRINCIPAL PLANE OF THE SUN 


In Part I, we also presented the observations of polarisation during 
twilight from two directions in the sky other than zenith, (1) at 30° from 
the zenith on the sunside 6 = 30°, ¢ = 0° and (2) at 30° from the zenith 
on the anti-sunside, @ = 30°, ¢ = 180°. These curves while showing the 
general features of the zenith polarisation curve, show some differences, 
In order to understand this, calculations have also been made of the polarisa- 
tion in these two directions. They are plotted in Fig. 11. 


The polarisation curves for the sunside show that the polarisation of 
primary scattered light (P) increases with Zp) and so also the polarisation of 
secondary scattered light (S). However, the secondary scattered radiation 
is polarised to a much smaller extent than the primary and hence the resultant 
polarisation P + S, increases with Zp) in the beginning when the primary 
scattered radiation is prominent and decreases afterwards when Z, = 94-5°, 
For Z, = 100°, the primary scattered radiation is negligible and P+§ 
practically coincides with S. On adding the night sky radiation, the P + § 
curve deviates from the S curve at 100°, and by 102° a sharp fall occurs, 


On the anti-sun side, the polarisation of P shows a fall with Z, and the 
rate of increase of S is much smaller than on the side of the sun. This causes 
a slight fall of polarisation in the beginning and a comparative flat region 
between 98° and 102°. 


VI. SUMMARY 


The formule for the intensity of the secondary scattered light from the 
sky as developed by Chapman and Hammad for a plane-parailel infinite 
atmosphere of defined optical thickness have been modified so as to deal 
with the twilight problem in which the curvature of the atmosphere is of 
fundamental importance. Calculations have been made of the intensity 
and polarisation of the light received from the zenith sky for different values 
of the sun’s zenith distance from 90° to 107° and for a radiation for which 
the optical thickness of the atmosphere is 0-2. For 0° to 5° solar depres- 
sion, the primary scattering is predominant; as the angle increases, the 
secondary scattered light becomes more and more important and when the 
sun’s depression is 8° to 12°, it is the principal contributing factor. Affter- 
wards, the night air-glow becomes comparable in intensity and when the 
solar depression is 17° to 18°, only the air-glow radiation remains. The 
calculated values of intensity and polarisation are compared with the observa- 
tional data collected at Mount Abu. 
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Fic. 11. Polarisation of the sky at points 30° from the zenith in the principal plane of the sun. 


.. Percentage polarisation of primary scattered light. 
. Percentage polarisation of secondary scattered light. 
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. Total effective percentage polarisation when night air-glow of constant 


- Observed sky polarisation at Mount Abu, 
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Effective percentage polarisation of primary plus secondary scattered light. 


intensity and 4% polarisation is added to primary plus secondary scat- 
tered light. 
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ABSTRACT 
Mineralogical compositions of five Indian clays have been deter- 
mined by the petrographic, X-ray, electron-microscopic and differential 
thermal analysis methods. The relative importance of the methods is 
discussed and the procedures are described. 


INTRODUCTION 


A stupy of the mineralogical composition of five samples of Indian clays, 
used in the ceramic and refractory industries, has been made. The import- 
ance of such a study lies in the fact that to a large extent the working proper- 
ties of the clays depend on their mineralogical composition, degree of crystal- 
linity and grain size. 

To determine the mineralogical composition of the clays, techniques 
complementary to one another, are generally required. As the grain size of 
most clay minerals is of the order of a few microns or less, the petrographic 
microscope does not have sufficient resolving power to distinguish the indi- 
vidual crystals. However, if the crystal size is coarse enough, as is the case 
with some dickites and kaolinites, the petrographic microscope can be used 
for their identification. On the other hand, it is possible to distinguish various 
types of clay minerals by X-ray diffraction procedures. The chief limitation 
of the X-ray method is that it is often difficult or impossible to detect impurities 
present in amounts less than 5% of the total material. It is here that the petro- 
graphic microscope has an advantage over the X-ray diffraction method. 


The electron microscope has been of great help in determining the mor- 
phology, degree of crystallinity and grain size, which aid in distinguishing the 
clay minerals. Differential thermal analysis has been of great value for the 
study of clay minerals, but here also, it is not always possible to detect small 
percentages of impurities present in the clays. In the present paper, the results 
are given of a correlated study using the above techniques. 





* This investigation was made by the author as part of his work on a Colombo Plan Fellow- 
ship at the Physical and Crystal Chemistry Section, Mineral Dressing and Process Metallurgy 
Division, Mines Branch, Department of Mines and Technical Surveys, Ottawa, Canada. 
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EQUIPMENT AND PROCEDURES 
Petrographic Study 


The clays were examined under a petrographic microscope by the powder 
method, using Shillaber’s oils of known refractive indices. The samples 
were examined as received and also after washing away the finer and lighter 
fractions. 


Electron Microscopic Study 


An R.C.A. electron microscope of the electromagnetic type was used in 
this study. The fine clay suspensions were mounted on “ formvar”’ films 
by standard procedures. The strength of the ‘“‘ formvar”’ solution used was 
0-1% in dichloro-ethane. The films were shadowed by chromium evapora- 
tion in vacuo. The samples were examined at an initial magnification of 
10,000. 


Differential Thermal Analysis 


Differential thermal analysis equipment has been described by many 
workers with some variation in details... A brief description of the experi- 
mental arrangement used in the present work is detailed below. Sketches of 
the palladium block specimen holder are shown in Fig. 1. 





DOR/ILL 2 HIOLES Z LIA. 3° DEEP N24/ DRILL DIA, 


— 
46 


Fic. 1. Palladium Specimen Hélder for D.T.A. 


The sample under study is placed in a cavity in the specimen holder and 
an inert reference material, calcined alumina, is placed in another cavity. 
Two connected thermocouples are placed one in each cavity and the leads 
of this differential couple are connected to a preamplifier and then to a Leeds 
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and Northrup Recorder. Another thermocouple is placed in a separate cavity 
close to the specimens. The leads of this couplé are connected to a Leeds 
and Northrup Program Controller. The horizontal tube furnace carrying 
the specimen holder is now heated at the rate of 12°C. per minute. The rate 
of temperature rise is controlled by the Program Controller. The inert mate- 
rial in the palladium block is stable throughout the heating range, whereas 
the clay sample undergoes a succession of endothermic or exothermic reac- 
tions during heating. This results in the setting up of differential e.m.f.s which 
are recorded on the Leeds and Northrup Recorder, appearing as peaks on 
either side of the centre line of a strip chart. However, as soon as each reac- 
tion is over the temperature of the specimen equalises with that of the inert 
material, and the differential thermocouple record returns to the centre of 
the strip chart. 


X-Ray Diffraction 


X-ray examination was carried out by the powder diffraction techniques 
using North American Philips Co. high angle goniometer. The technique 
is described in detail by Klug and Alexander.? The salient features are as 
follows : 


An X-ray beam impinges on the smooth surface of a finely pulverised 
specimen mounted in a rotating sample holder. The diffracted X-ray beam 
is picked up by a Geiger-Muller counter, scanning at a rate of 2 9 = 1° per 
minute, amplified and recorded on a moving chart from which @ values can 
be obtained and hence the corresponding d-spacings. 


RESULTS 

X-Ray Diffraction 

The results of X-ray examination of these five clays are given in Table I. 
Filtered Cu K radiation was used throughout. The d-spacings in Angstrom 
units are given for each clay. By comparison of the values with those in the 
A.S.T.M. Card Index of Powder Diffraction Patterns and also the available 
lieterature,* each set of the lines was attributed to a specific mineral consti- 
tuent. The intensities of the kaolinite and mica lines showed some variations 
from those given by some of the previous workers. However, this could be 
expected, as the clays and micas do tend to exhibit preferred orientation. The 
intensities are estimated approximately from the magnitude of the peaks. 
Quartz, even when present in small amounts, gives much stronger peaks than 
kaolinite as, for instance, in the Rajhara and Jabbalpur clays. In such cases 
the strongest peak of the kaolinite is taken as 100, with which other peaks of 
lower intensities are compared, while the quartz lines of higher intensity are 
denoted as “‘ very strong’. The X-ray diffraction patterns of these clays at, 
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reproduced in Fig. 3. These patterns were recorded for reproduction pur- 
poses with the recorder chart moving at the rate of 9” per hour. These records 
show a slightly more irregular background than the originals, from which the 
d-values were calculated, which were recorded at 36” per hour. 


The differential thermal curves of these clays are given in Fig. 2. All 
these clays show broad endothermic peaks between 586-610° C. and a com- 
paratively sharp exothermic peak between 975 and 990° C., characteristic of 
kaolinite. The endothermic peak is thought to be due to the loss of combined 
water and the exothermic peak due to the crystallisation of alpha alumina 
or mullite, or both A small exothermic peak at 450°C. in the Rajhara 
clay appears to be due to the oxidation of some organic matter. During 
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Fic. 2. Differential Thermal Curves. 
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cooling all the clays showed a small exothermic peak at 573°C. which is 
attributed to the 8 to a inversion of quartz. 
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OEGREES 2806 
Fic, 3. X-Ray Diffraction Patterns of Indian Clays. 


The electron micrographs of these clays are shown in Figs. 4-8. The 
Travancore-China clay is extremely well crystallised, and shows equant or 
elongated crystals of six-sided cross-section. The Jabbalpur clay also shows 
many such well-formed crystals besides the coarsely formed irregular plates. 
The Orissa clay shows a few subhedral crystals with the irregularly formed 
plates predominating. The Rajhara clay is mostly composed of irregularly 
formed coarse plates with some smaller subhedral to anhedral crystals. The 
shape of the crystals in all these clays is suggestive of kaolinite with a possi- 
bility of the coarser plates being dickite. A recent publication’ shows that 
kaolinite could occur in tubular forms. However, such forms were not seen 
in the present study. The Chaibasa clay is very coarsely crystalline, showing 
large irregular plates. Some of these could possibly be mica. 


Under the petrographic microscope the clay mineral, in most cases, 
appears as aggregates of minute plates and needles, with an average refractive 
index of 1-562, and with moderate to weak birefringence. In the Chaibasa 
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clay, some of these small plates, under high magnification, appear to be vermi- 
cular. The impurities present in these clays are as follows: 


Travancore-China clay A small amount of quartz with traces of rutile, 
zircon, apatite, hydrated iron oxide and musco- 
vite. 

Chaibasa clay Appreciable amounts of muscovite and kyanite, 
small amounts of serpentine and quartz, and a 
trace of hematite. 

Orissa clay Small amounts of quartz and muscovite, with traces 
of rutile, limonite, apatite, dolomite, zircon, 
sericite and orthoclase. 

Jabbalpur clay Small amounts of quartz and sericite (muscovite) 
and talc, with traces of rutile, apatite, calcite and 
gibbsite. 

Rajhara clay A considerable amount of quartz with a small 
amount of talc, traces of apatite, calcite, and hy- 
drated iron oxides. 


SUMMARY 


The principal constituent in all these clays is kaolinite. It is fine-grained 
and very well crystallised in the Travancore-China clay. In Jabbalpur, Orissa 
and Rajhara clays, the kaolinite occurs mostly as coarse irregular plates with 
some fine euhedral to subhedral crystals. The Chaibasa clay shows the 
coarsest plates, some of which seem to be muscovite. Amongst the impuri- 
ties, quartz occurs in all these clays in varying amounts, the Rajhara clay 
containing the largest amount. Muscovite occurs in a significant amount in 
the Chaibasa clay, and in small amounts in Orissa, Jabbalpur and Travancore- 
China clays. Kyanite occurs in an appreciable amount in Chaibasa clay, 
along with a small amount of serpentine. Small amounts of talc occur in 
the Rajhara and Jabbalpur clays. The Rajhara clay contains some organic 
matter. Other impurities occurring in very small amounts or traces are rutile, 
hematite, apatite, zircon, orthoclase, calcite, dolomite, gibbsite and hydrated 
iron oxides. 
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Fic. 6. Orissa Clay, xX 20,000 
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THE ROLE OF SOLVENT IN CHEMICAL 
REACTIONS 


Part II.* The Reaction between Aniline and Methyl Iodide 
in Different Solvents 


By V. S. PADMANABHAN{ 
(Contribution from the Chemistry Department, Madras Christian College, Tambaram) 
Received March 19, 1956 


(Communicated by Dr. S. V. Anantakrishnan, F.A.sc.) 


IN the previous part! were reported some observations on the influence of 
solvents on the reaction between pyridine and methyl iodide. The present 
studies deal with the alkylation of a primary amine. While the first reaction 
was comparatively simple, here there are two complications—consecutive 
reactions of progressive alkylation, replacing the hydrogen atom of the 
primary amine and then formation of the quaternary ammonium compound. 
While the two substitution reactions may involve comparable activation 
energies, the formation of the quaternary ammoium compound may involve 
a lower activation energy. Cox? has studied the reaction between aniline and 
w-bromo acetophenone in several solvents while Hawkins* has studied the 
reaction between pyridine and allyl bromide in different solvents. Neither 
worker has been able to notice any specific correlation between any of the 
physical properties of the solvent and the reaction rate. The reactions had 
also been studied by these workers under conditions in which heterogeneity 
could not be excluded. Subsequent work also indicates the use of benzyl 
bromide* and allyl chloride® in the reaction with pyridine and it has been 
observed that the influence of dielectric constant is sometimes in opposite 
directions when the different halides are compared. It was felt that the use 
of methyl halides in the first instance will simplify the situation and serve as 
a starting point in the analysis of solvent influences. 


The reaction between aniline and methyl iodide proceeded under homo- 
geneous conditions in both acetone and ethyl alcohol but in nitrobenzene and 
carbon tetrachloride solutions even in very dilute solutions separation of a 
solid product was evident quite early in the course of the reaction. The 


* Based on the Thesis approved for Ph.D. Degree of Madras U-:iversity. 
~ Now at the Central Leather Research Institute, Madras-20. 
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analysis of both precipitate and solution from time to time for nitrogen as 
well as iodide contents showed clearly that the solid separating was primarily 
aniline hydriodide while the solution contained the products of progressive 
alkylation, depending on the duration of the reaction. The reaction apparently 
did not go to completion, the amount of methyl iodide used up being only 
about 87% at about 35° C. even after eight months. 


Irrespective of the solvent used, the initial rate, whether one takes the 
value obtained by graphical extrapolation of experimental data or actual 
calculated value for the first 10% of the reaction, differed to varying extents, 
from the mean value of the rate constant. The reaction was found to be best 
considered as a second order one even though there was a pronounced drift 
in the value. For purpose of discussion, the values are collected together in 
Tables I, I, III and the drift in the values of the rate constants illustrated 
in Figs. 1 and 2 below: 


TABLE I 


The reaction in nitrobenzene 





[Aniline] : Tempe- Initial value of Mean value of 
[methyl rature bimolecular rate bimolecular rate 
iodide] a ot constant x 104 constant x 104 





1-50 50 2°33 -80 
1-53 60 4-29 69 
1-15 40 1-62 “Fa 
50 2°29 63 
45 1-62 “88 
50 2°31 a? : 
50 2°45 
55 3-46 
60 5-3 
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TABLE II 


The reaction in acetone 





[Aniline] : Tempe- Initial value of Mean value of 
[methyl rature bimolecular rate bimolecular rate 
iodide] °C. constant x 104 constant 104 





1-44 40 2-70 3-11 
1-28 5-74 6-03 
0-79 1-87 1-64 
0-88 2°37 2°41 
0-88 : 3-30 3-35 
0-91 4-8] 4-50 
0-86 4-22 4-57 
0-87 6-11 6°18 





TABLE III 





Solvent Temperature® C. % Reaction k, x 104 


Nitrobenzene a 40 8-49 1-618 
27°57 1-646 
43-03 1-706 
55°76 1-738 
63-95 1-793 





95% Ethyl alcohol .. 28 -39 16-82 
48-12 17-55 
62-10 17-51 


Acetone .. si 13-11 4°81 
29-38 4°75 
39-62 4-55 
47-37 4-49 
57-10 4°35 





N.B.:—In all the tables, concentrations are approx. M/10, all rate constants in moles. 
litre and the unit of time is the second. 
ja 


2. In Table III, k,-First order rate const. calcd. with respect to methyl iodide. 
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Fic. 2. Reaction in Acetone 
Grant and Hinshelwood® in their study of benzoylation of aniline found 
that the best constants were obtained when the proportion of amine to acyl 
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halide was 2:1, the second molecule of the amine serving to prevent the 
reversal of the reaction. If similar conditions were to hold good in the alky- 
lation studied here, with the consumption of methyl iodide as the measured 
quantity, one might expect that the rate expression of the type: 


dx/dt = k,(a — x) (b — 2x) 
will apply but calculations showed that no constant could be obtained using 


the integrated form of this equation. The formation of aniline hydrochloride 
as well as the drift in the rate constant call for a different explanation. 


The following stages can be visualised for the reaction: 
C;H;NH, + Mel ——>(C,H;NH.Me)* + I- (i) 
(CsH;NH,Me)* +>C.H;NHMe + H+ (ii) 
(CsH;NH,Me)t > C;H;NH, + Met (iii) 


Ionisation of methyl iodide is unlikely to be first stage in reaction since the 
rate is found to be quite dependent on the initial concentration. The first 
reaction is clearly the rate-determining step and the general trend of rate con- 
stants is that of a second order. (ii) and (iii) are reversible reactions and in 
the presence of an adequate quantity of aniline, reaction (iii) will be negligible. 
The reversibility could be readily demonstrated by starting with an alkyl 
aniline hydrochloride. After contact with a solution of aniline in the differ- 
ent solvents, the final product analysed to indicate the interchange, the solid 
being essentially aniline hydrochloride. 


TABLE IV 





Bimolecular rate Arrhenius 
Dielectric constants x 104 parameters 
Solvent constant 
at20°C. 40° 45° 50° 55° 60° E logyoPZ 
€. €. Cc. ¢€. c. K. Cals. 








Acetone + 2:17 3:09 4-06 5-88 


Ethyl alcohol 
(95%) ~ oa .. Pe . Be 15-37 7°35 


Nitrobenzene 35-79 a 1-60 2:44 3:46 5-67 16-83 7-80 





N.B.—The rate constants in acetune and nitrobenzene have been obtained from experiments 
in which the ratio Aniline: Methyl iodide was in the region 0-85-0-90. 


Since the reaction appeared to be sensitive to even small changes in the 
concentrations, the order of the reaction could be computed only by van't 





The Role of Solvent in Chemical Reactions—II 373 


Hoff’s differential method and these indicate that the reaction is of the first 
order with reference to aniline but the reaction order with respect to methyl 
iodide is only half. Since the measured quantity is only the change in methyl 
iodide concentrations and part of the iodide is fixed as aniline hydriodide, 
further work is needed to indicate the path of the reaction in full. With this 
limitation, however, we can analyse the influence of the solvent. 


For reasons already indicated, comparison could be made using the 
initial value of the bimolecular rate constant (obtained by extrapolating 
k,-% reaction graph to 0% reaction), in the different’ solvents at different 
temperatures. Some justification for this procedure’is the observation that 
the log k — 1/T graph showed’ very good linear relationship (Fig. 3). 


It is clear that while there is a corresponding’change in both Arrhenius’ 
parameters similar to that indicated by Fairclough and Hinshelwood,’ there 
does not appear to be any correlation with dielectric constant, of the type 
expected for a dipole-dipole reaction. Solvation of the transition state, 
removal of the proton as a solvated ion in the second stage of the reaction, 
changes in the dielectric constant with temperature (which have not been 
taken note of in the present calculations), solubility of the different products 
—all are relevant factors to be taken into account and further information is 
clearly needed. As a first step, solvent mixtures might be expected to give 
useful data and these are considered/in the next part. The complex nature 
of specific solvent influences are shown clearly by the rate constant, the reac- 
tion being fastest with 95% ethyl, alcohol even though the activation energy 
in this solvent is much higher than in acetone. 


If the solvent is to take part in proton removal, the alternative question 
arises: Is there any unimolecular change with ionisation of methyl iodide 
as the rate determining step? This step cannot be ruled out on account of 
the high dielectric constant of the media. Hughes, Ingold and Shapiro® have 
shown, for instance, that where, as in the case of Isopropyl halides, both 
§,1 and S,2 mechanisms are possible, the rate constant may show a continu- 
ous rise. This feature is present in the present reaction at the lower tempera- 
tures but at the highest temperatures studied and probably at the later stage 
in the reaction the trend is reversed and there is a sign of decrease. At the 
present stage, it is not possible to isolate the two types of reaction and with 
a primary amine further alkylations which also consume the alkyl halides, 
are unavoidable so that the consumption of methyl iodide appears to be more 
than can be accounted for by the first stage of the alkylation. In either case, 
it is clear that the introduction of an alkyl group is analogous to the formation 
of ‘ Onium compounds ’, 
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L SOLVENT NITROBENZENE 
Il SOLVENT ACETONE 





2 





2.98 3.00 a1 


+ x10? 


Fic. 3 
EXPERIMENTAL 


Aniline, redistilled, B.D.H., was distilled once, kept over Analar potas- 
sium hydroxide for a week and then redistilled with the minimum exposure 
to light. Even diffuse daylight was found to cause the deepening in the colour 
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of the distillate. Rejecting large head and tail fractions, only the middle 
fraction was collected, B.P. 183° C./760 mm. 
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Methyl iodide was purified as in the previous part, 


Nitrobenzene, pure redistilled, B.D.H., was shaken small with 
small portions of sodium carbonate solution till the aqueous portion was 
colourless, then washed repeatedly with distilled water, left over anhydrous 
calcium chloride in the dark for a week and then fractionated, only the middle 
fraction being used, B.P. 210° C./760mm. Moelwyn-Hughes® considered 
that this drying agent was not adequate but it was found that by using long 
periods of drying and rejecting large head and tail fractions, contamination 
by moisture is negligible. The amine hydriodides formed in the reaction 
are found to be appreciably soluble even in the presence of minute amounts 
of water in the solvent and a test for the absence of water in the solvents is 
the ready precipitation of the salts even in the most dilute solutions. Further, 
it was observed that freshly distilled nitrobenzene is very - in colour 
which deepens with the absorption of moisture. 


Acetone of Laboratory Reagent quality was purified as in the previous 
part.! 


Carbon tetrachloride of redistilled grade was left over mercury, shaken 
at frequent intervals until there was no further formation of mercuric sulphide, 
then washed thoroughly with sodium bicarbonate and then with distilled 
water, and finally dried over anhydrous calcium chloride and fractionated 
B.P. 75-6° C./760 mm. 


Preliminary experiments.—Aniline and methyl iodide reacted vigorously 
with considerable evolution of heat. In nitrobenzene, the reaction was still 
exothermic but much slower. Leafy white crystals with no sharp melting 
are generally obtained. Estimation of the iodide content gave a value 
nearer to aniline hydriodide than to the alkylated product. (In the 
course of these studies, it has been found that estimations using Fajans’ 
method with eosin as the indicator gave more consistent results.) 
The results show clearly that the product separating out was essentially 
aniline hydriodide with about 3% of methyl aniline hydriodide. Nitrogen 
estimations by microkjeldahl methods gave identical results. Relative 
proportions did not alter the nature of the precipitated product. The 
examination of the mother liquor indicated the presence of alkylated products 
at all stages of methylation. The reaction in carbon tetrachloride was slower 
and the precipitate obtained had as before no sharp melting point. Esti- 
mation of iodide content indicated that invariably it was of mixed composi- 
tion. Similar results were obtained using chloroform and acetic acid as 
solvents, the reaction being extremely slow in the last named solvent, 


Typical results of analyses are presented in Table V, 
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TABLE V 





Nitrobenzene Carbon tetra- 
solutions chloride solutions . 





Precipi- Mother 
tate liquor 





Experiment No. 45 11 


50 
(Left over for 8 
months) 


Concentration of aniline... M/10 M/5 M/5 M/10 M/5 


Concentration of methyl 
iodide Nt .. M/10 M/10 M/5 M/10 M/5 


Per cent. nitrogen . GR 6-305 
Per cent. iodine - ee 48 -67 56°15 56:08 & 56-74 





Aniline hydriodide requires nitrogen 6-339% and iodine 57-46%. 


Kinetic studies—The reactions were studied in a large thermostat, elec- 
trically operated and maintained at the requisite temperature with a tolerance 
of + 0-005 at 40°, + 0-01 at 50° and + 0-015 at 60° C. The reactant 
solutions were prepared by weighing both solvent and solute and thoroughly 
mixing the two in well-stoppered bottles. Measured equal volumes of 
these were pipetted into bulbs of borosilicate glass and sealed. A number 
of these sealed tubes were kept in the thermostat, and from time to time 
taken out and analysed. 


The contents of the bulb were transferred to a separating funnel, rinsing 
the bulb with ice-cold distilled water, shaken with benzene sufficient to form 
an upper layer dissolving the reaction mixture. The aqueous layer was tapped 
off and the benzene layer extracted four times with ice-cold water. All the 
aqueous portions were combined and the iodide content estimated by 
Fajans’ method using eosin as indicator. 


Where the solvent was miscible with water, the contents of the bulb 
were taken up with ice-cold water to which sufficient acetic acid was added 
to chill the reaction and the iodide estimated with silver nitrate solution. 


Typical runs in the different solvents are presented in Tables VI to IX, 
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TABLE VI 
Reaction in nitrobenzene 





Temperature 40° C. 50° C. 60° C. 





Aniline concentration 0-09434 mol./I. 0-03466 mol./l. 0-1074 mol./l. 


Methyl iodide con- 
centration .. 0-08414 ,, 0:01839 0-07004 sé, 


Strength of AgNO, 
solution ar 0-0255 N 0-01048 N 0-02462 N 





Time V Time Time V 





5880 2°75 19780 . 3600 4-15 
18876 7°55 43890 ° 7200 7-95 
31600 11-15 49855 , 10800 10-45 
74830 18-85 59800 ° 14400 12-65 

103500 21-85 77650 ° 18000 14-98 
132400 24-20 95940 > 21600 16-75 





TABLE VII 
Reaction at 50° C. in nitrobenzene for elucidating reaction order 





Concentration 0-09674 mol./l. 0-09674 mol./l. 0-09731 mol./l. 0-1390 mol./l. 
of aniline 


Concentration 0-09491 mol./l. 0-13060 mol./l. 0-1098 mol./l. 0-1098 mol./l. 
of methyl io- 
dide 


Strength of 0-02514 N in all the four experiments 
AgNO, solu- 
tion 
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TABLE VIII 
Reaction in 95% ethyl alcohol 





Cee as 


Temperature °C. 40 50 60 





Concentration of ani- 0-09408 mol./I. 0-08732 mol./I. 0-1075 mol./I. 
line 

Concentration of 0:08553_—,, 0:08839 ss, 0-08543_—s,, 
methyl iodide 


Strength of AgNO, 0-02550 N. 0-02504 N 0-02462 N. 
solution 





V Time 





7:90 2100 
12-60 4500 
16-10 7500 
18-60 9300 
20-30 11700 
21-6 14100 
24-3 16500 





TABLE IX 
Reaction in acetone 





Temperature °C. 40 50 





Concentration of aniline Pe 0-08935 mol./l. 0-08650 mol./lI. 
Concentration of methyl iodide 0-08946__—sé,, 0-08688 __,, 
Strength of AgNO; solution 0-02550 N. 0-02496 N. 





Time Time 





9218 7825 
16460 ; 14400 
23570 25200 
32700 39600 
47130 61200 
76320 





Time in Seconds; V, Titre value in ml. 
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Examples of the general trend with other concentrations are illustrated in 
Figs. 5 and 6. 


TATOcC WALIEC tee Mee 


a 
z 
= 
w 
2 
J 
< 
> 
wi 
io 4 
= 
K 


TEMPERATURE 50°C 
A,,Az METHYL 10DIDE CONCENTRATION SAME, 


ANILINE A,:A2:.1°4:5 
GB). By, ANILINE CONCENTRATION SAME, 


METHYL IODIDE B,.B2:. 1°15 





200 4.00 
TIME IN MINUTES 





Fic. 5. Solvent Nitrobenzene 
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TITRE VALVE IN ML 








200 
TIME IN MINUTES 


Fic. 6. Reaction in acetone 


No kinetic studies were made in pure carbon tetrachloride as it was not 
possible even in the most dilute solutions to have homogeneous conditions. 
Similarly no kinetic studies were made in glacial acetic acid as the extent of 
reaction was small even after weeks. 


SUMMARY 


The reaction between aniline and methyl iodide has been studied at four 
different temperatures in 95% ethyl alcohol, acetone, nitrobenzene and carbon 
tetrachloride. The complications arising from progressive alkylations as 
well as from the crystallisation have been minimised under the experimental 
conditions in the initial stages. The observed changes in the rate constants 
and in the Arrhenius parameters are discussed. 
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INTRODUCTION 


Jessop! has mentioned that cadmium gives rise to polarographic wayes in 
presence of triethanolamine, but states that the top of the steps are rounded. 
Work conducted in this laboratory, however, has indicated that very good 
waves are obtained in presence of ethanolamines under certain conditions. 
In the present communication a detailed account is given of the polarographic 
behaviour of cadmium in presence of ethanolamines. 


EXPERIMENTAL 


The procedure employed in this work is essentially the same as that 
described previously.2 The values of ‘m’ for the two dropping mercury 
electrodes used are 1-227 and 1-150 mg. per second respectively, the drop- 
time ‘t’ being five seconds in normal potassium chloride. To facilitate 
comparison all the wave-heights recorded have been reduced to the corres- 
ponding values of the capillary with m= 1-150mg. per second. 


RESULTS 
1. Effect of pH 


The results are summarised in Table I. In solutions containing mono- 
and diethanolamines at pH below 11, the polarographic wave seems to 
consist of two to three waves, the half-wave potentials of which are very 
close to each other. Beyond pH 11, these wavelets disappear and only 
one reversible wave appears. In triethanolamine solutions, between pH 9-8 
and 6-5 two waves are noticed (Fig. 1). In more alkaline solutions, however, 
the second wave appears in the form of an elongated bend before the 
diffusion current plateau (Fig. 2). As the pH of the solution is increased, 
the half-wave potential increases slowly till a pH of about 9-8 is reached 
and thereafter it increases very rapidly. The value of E3,,— E,,, varies 
between 0-031 and 0-049 V at different pH values indicating both reversible 
and irreversible electrode processes. 


AS 383 
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2. Effect of variation of the composition of the base electrolyte 


The effect of addition of sodium hydroxide, sodium carbonate and 
ammonium chloride-ammonium hydroxide as the reagents to alter the pH has 
been investigated and the results obtained are indicated below. It may be 
pointed out that the half-wave potentials are not altered appreciably in 
solutions of pH less than 9. 


TABLE I 


Effect of pH on the Polarographic Behaviour of Cadmium in Ethanolamines 





Half-wave Half-wave potential 
potential E,,4,—Ej,/4 vs. S.C.E. 
Amine used vs. S.C.E. Ist wave 
Ist wave (Volt) 2nd wave 3rd wave 
(Volt) (Volt) (Volt) 








Monoethanolamine .. , —0-57 - —0-60 
—0°63 ‘5 —0:70 
—0°65 oe —0-70 
—0-73 as 
—0-78 


Diethanolamine es : —0:57 
—0-61 
—0-65 
—0-75 
—0-80 


Triethanolamine a . —0-597 0-031 
—0-609 0-031 
—0-613 0-033 
—0-634 0-039 
—0-645 0-043 
—0-720 0-040 
—0-771 0-03] 
—0-870 0-049 





(i) Sodium hydroxide—A study of the polarograms obtained under 
different conditions has shown that in solutions containing triethanolamine 
the elongated bend (Figs. 2 and 3) before the diffusion current plateau is 
more marked in solutions with low amine/sodium hydroxide ratio, while 
in solutions containing mono- and diethanolamines perfectly sharp 
bends occur. The analysis of the waves obtained under varying conditions 
are presented in Table Il. The half-wave potentials and the E3,4— Ey, 
values of the waves for the triethanolamine complex are reported after 
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applying correction for the bend and are hence only approximate. From 
Table II it is clear that at a constant concentration of sodium hydroxide 
the half-wave potentials and E;,,—E,,, values increase with an increase in 
the concentration of triethanolamine. It should be noted that the value 
of Es4— Exa is about 0-035 V. at low concentrations of triethanolamine 
corresponding practically to a reversible 2-electron reduction step. In the 


TABLE II 


Effect of Sodium Hydroxide on the Polarographic Behaviour of Cadmium 
in Ethanolamines 





Concentration of Half-wave Concentra- 
potential Height tion of 
Amine used Amine Sodium vs. E3;4—E,,;4 of the cadmium 
(Molar) hydroxide _ S.C.E. (Volt) wave S/30 (Milli- 
(Molar) (Volt) molar) 








—0-753 0-027 _ slight pptn. 
—0-777 0-032 9 
—0:794 0-029 47-5 
—0-815 0-025 _ slight pptn. 
—0-836 0-032 se 
—0-849 0-032 44-0 
—0:750 0-028 © slight pptn. 
—0-773 0-029 45-7 
—0-787 0-027 43-6 
—0-324 0-030 42-2 
—0-846 0-033 41-8 
—0-857 0-030 38-8 
—0-747 0-042 sia 
—0-771 0-031 49-1 
—0-825 0-051 46-2 
—0-857 0-124 43-8 
—0-885 0-099 41-3 
(Small bend) 
0:10 —0:936 06-114 37-0 
(Bend has almost disappeared) 
0-50 —0-837 0-036 Me 
0-50 —0-859 0-047 48-4 
(Slight elongated bend) 

0-50 —1-:010 0-102 35-2 

. —0-860 0-036 * 
—0-870 0-049 46-2 
—0-929 0-087 43-3 
—0-976 0-099 40-2 
—0-995 0-102 37°8 
—1-036 0-093 31-0 


Monoethanol- 
amine 


Diethanolamine 


Triethanolamine 


SOSOTOOTOOTSOTOO 
SSSSSASSASSASSAS 
eeesers eters ees 
SSsssSSSsSssSSSsssq 
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case of solutions containing mono- and diethanolamines, the half-waye 
potentials increase with the increase in the concentration of the amine while 
the value of E;,,— E,,,4 is of the order of 0-028 V. under all conditions ind. 
cating a reversible 2-electron process. The wave-heights given in Table II 
indicate a decrease in the wave-height with an increase in the number of 
ethanol groups in the amine. 


(ii) Sodium carbonate.—The results presented in Table III below indi. 
cate that the general behaviour in presence of sodium carbonate is the same 
as in the case of sodium hydroxide. 


TABLE III 


Effect of Sodium Carbonate on the Polarographic Behaviour of Cadmium 
in Ethanolamines 


Concentration of cadmium = 1-089 millimolar 





Concentration of 





Half-wave 
Sodium potential E;,,—E,/, 
carbo- vs.S.C.E. (Volt) 
nate (Volt) 
(Molar) 


Amine used Height of the 


Amine wave 


(Molar) 





—0-739 
—0-764 
—0-772 
—0-767 
—0-787 
—0-797 
—0-710 


0-029 
0-032 
0-032 


Monoethanol- ; 
0-029 : 


0 
amine 4 
4 
0-029 
0-023 
0-028 


Ss8S88sss5 


Diethanolamine 


ptn.) 
6 


—0-731 0-033 


Triethanolamine 


Se S32 T9 Ort? ree 
SsSassa ssassas 


S385 S S38 


-OOo°O 


s8ssss 


Cae Ceo OorrrsSe or rr ooo 
So SO 


ssss 


—0-738 
—0-746 
—0-759 
—0-770 
—0-707 


—0-739 
—0-767 
~ 0-915} 
—0-806) 
—0-951f 
—0-759 
—0-789 
—0-822 
—0-859 


0-031 
0-030 
0-030 
0-031 
0-037 


0: 
rere 
0-109 
rite 
0-100 
0-055 
0-058 
0-063 
0-800 


(Slight 
45: 
42: 


3 
8 
3 
8 
Pp 
6 
3 


40- 
36°6 
32°7 
38+1 


(Slight precipitation) 
44-7 


41-3 
37°6 


45-2 
35-3 
31-9 
28-6 
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(iii) Ammonium chloride (1 M)—Ammonium hydroxide (1 M).—The 
waves are very well defined. The half-wave potentials, the values of 
Ey«— Exya, and the heights of the waves are given in Table IV. It is clear 
that the half-wave potentials increase with an increase in the concentration 
of the amine. The value of E;,,— E,,4 increases with an increase in the 
concentration of triethanolamine but in the case of mono- and diethanol- 
amines the values are round about 0-029 V. indicating a reversible two- 
electron reduction process. The height of the wave diminishes very rapidly 
with increasing concentration of triethanolamine whereas in the case of mono- 
and diethanolamines the fall is not so great. 


TABLE IV 


Effect of Ammonium Chloride-Ammonium Hydroxide on the Polarographic 
Behaviour of Cadmium in Ethanolamines 





Concen- Half-wave Height Concentra- 
Amine used tration of potential E;,4,—E,,4 of the tion of 
amine vs. S.C.E. (Volt) wave cadmium 

(Molar) (Volt) at S/30 (Millimolar) 





Monoethanolamine —0-784 0-027 
—0-791 0-030 


—0-815 0-028 


“OO 
Ss8 


—0-794 0-029 
—0-795 0-029 
—0-807 0-026 


Diethanolamine 


-_OOS 


—0-757 0-029 
—0-762 0-033 
—0-772 0-043 57°5 
—0-786 0-041 48-4 
—0-798 0-041 45-4 


Triethanolamine 


eee 
SSSs8 Ss8 





3. Effect of concentration of cadmium on the diffusion current constant in 
presence of sodium hydroxide, sodium carbonate and ammonium 
chloride (1 M)-ammonium hydroxide (1 M). 


The results obtained are presented in Table V. In solutions containing 
sodium hydroxide or carbonate, 0-5 or 1-0 M amine is employed due to low 
solubility of the complex. It has been noticed that the solubility of the 
monoethanolamine complex in sodium carbonate is greater than in sodium 
hydroxide solutions in contrast to the behaviour of the diethanolamine 
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TABLE V 


Effect of Concentration of Cadmium on the Diffusion Current Constant* 





Diffusion 
Conc. of current 
cadmium constant 
(Millimolar) ig 
Cm 2/3 p 16 


Diffusion 
current 
constant 


Concn. of 
cadmium 
(Millimolar) 





(1) B.S.: 1-00 M Dien+-1-00 M NaOH 
X = 1-436 at — 1-0 V. 
1-089 2-32 
4-356 2°32 


(2) B.S.: 0-50 M Trien+0:10 M NaOH 
X = 1-417 at — 1-3 V. 
1-089 2°51 
5-445 2°49 

(3) B.S.: 0-50 M Trien+-1-00 M NaOH 
X = 1-417 at — 1-3 V. 
1-089 2°31 
5-445 2°28 


(4) B.S.: 1-00 M Moen+1-00 
Na.CO, 
X = 1-435 at — 1-0 V. 
1-089 2-35 
5-445 2°33 
10-89 2-32 
(5) B.S.: 1-00 M Moen-+-0-10 M 
Na,CO, 
X = 1:443 at — 1-0 V. 
1-089 3-05 
10-89 2:71 
(6) B.S.: 1-00 M Trien+-0-10 M 
Na,CO, 
X = 1-426 at —1-1V. 
1-089 2-29 
5-445 2-24 


(7) B.S.: 1-00 M Trien+-1-00 M 
Na,COgs 
X = 1-426 at — 1-1 V. 
1-089 1-74 
5:445 1-74 
(8) B.S.: 0-10 M Moen+ Ammonium 
chloride (1M)+ Ammonium hydro- 
xide (1 M) 
X = 1-447 at —1-0V. 
1-089 3-96 
5-445 3-94 
(9) B.S.: 0-10 M Dien+ Ammonium 
chloride (1 M)+- Ammonium hydro- 
xide (1 M) 
X = 1-447 at —1-0V. 
1-089 3-90 
5-445 3-88 
(10) B.S.: 0-02 M_ Trien+ Ammonium 
chloride (1 M)+- Ammonium hydro- 
xide (1 M) 
X = 1-447 at —1-0V. 
1-034 3-96 
5-170 3-93 
(11) B.S.: 0:-60M_ Trien-+- Ammonium 
chloride (1 M)+ Ammonium hydro- 
xide (1M) 
X = 1-431 at —1-0V. 
1-034 3-05 
5-170 3-03 





* Abbreviations :—B.S.-Base solution ; X-m*/*r/6 in mg.?sec.—/2 ; Moen-monoethanol- 
amine ; Dien-diethanolamine ; Trien-triethanolamine. 


complex whose solubility is greater in sodium hydroxide solutions. 
bility of the monoethanolamine complex in sodium hydroxide and that of 
diethanolamine complex in sodium carbonate is far too low to be of 
practical value. The results indicate that there is good proportionality 


Solu- 
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between diffusion current and the concentration of cadmium. The con- 
stancy of the diffusion current constant in presence of ammonium chloride- 
ammonium hydroxide is quite good, indicating that the presence of ethanol- 
amines has no adverse effect. 


DISCUSSION 


1. Splitting of the polarographic wave in triethanolamine 


In solutions containing 0-1 M triethanolamine, it has been noticed that 
a doublet first appears at a pH 6-5, the half-wave potentials of the first and 
the second waves being ca— 0-61 V. and ca— 0-80 V. vs. S.C.E. respectively. 
As the pH is increased, the half-wave potential of the first wave increases 
(Table I) while that of the second remains practically constant. The doublet 
almost disappears as the pH is increased beyond 9-8. These observations 
can be explained as follows:—The single wave that appears at pH 5-35 is 
that of aquo-cadmium ion. As the pH is increased to 9-8, the first wave 
practically corresponds to that of aquo-cadmium ion while the second wave 
appearing at abou t— 0-75 V. vs. S.C.E. is perhaps due to triethanolamine- 
cadmium complex. When the pH is increased beyond 9-8, the half-wave 
potential of the first wave increases rapidly and when it approaches a value 
of ca — 0-75 V. vs. S.C.E. only a single wave with very little overvoltage is 
obtained. The wave-splitting is not due to two reduction states as the reduction 
of cadmium at the dropping mercury electrode involves only the straight 
two-electron reduction step, but due to the reduction of two species of com- 
plex ions in sluggish equilibrium. Such results have also been noticed by 
Lingane* while studying the polarographic behaviour of iron in tartrate 
media between pH 6 and 8. 


When the base solutions contain 0-1 M triethanolamine and sodium 
hydroxide (Table Il) or carbonate (Table III) the reduction process is very 
nearly reversible. As the concentration of triethanolamine is increased, 
the bend before the diffusion current plateau becomes elongated (at the same 
time E,,, — E,,, value increases) indicating the formation of another species 
(containing more of the triethanolamine) in sluggish equilibrium. 


2. Splitting of the polarographic wave in mono- and diethanolamines 


In contrast to the behaviour in triethanolamine, the splitting of the 
polarographic wave is practically absent in base solutions containing mono- 
and diethanolamines except below pH 11. Below pH 11 two or three poorly 
separated waves (Table I) are seen, the lowest wave corresponding to aquo- 
cadmium ion while the next two are due to the complexes containing (i) cad- 
mium ion + amine, (ii) cadmium ion + amine + hydroxyl. The formation 
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of the last species starts at pH 9-45. When the pH is increased beyond 11, 
only one wave is obtained showing the formation of complex (ii). 


3. Relationship between the height of the polarographic wave and the number 
of ethanol groups present in the complex-forming amine 


It has been noticed (Tables II, III, [V) that the wave-height for a given 
concentration of cadmium is highest in the case of monoethanolamine and 
falls off with an increase in the number of ethanol groups in the complex- 
forming amine. This indicates that the diffusion coefficient diminishes 
with an increase in the number of ethanol groups in the amine. 


4. Complex formation in presence of sodium hydroxide 


In base solutions containing mono- and diethanolamine, cadmium gives 
a reversible two-electron reduction step. Applying Lingane’s? equation to the 
data presented in Table II we get that the complexes have the formule 
[Cd (Moen), (OH),], [Cd (Moen); (OH),] and [Cd (Dien), (OH),]. In this 
connection it has to be pointed out that the shift in the half-wave potential for 
a given concentration range of the amine is always smaller in solutions con- 
taining higher concentrations of sodium hydroxide. This observation may 
be explained on the basis of the co-ordination taking place through the 
hydroxyl of the ethanolamine. When there is very high concentration of 
the hydroxyl group the half-wave potential will have been shifted to a suffi- 
ciently negative potential which diminishes the effect of co-ordination through 
the hydroxyl of the ethanolamine. 


The behaviour in triethanolamine solutions is completely different from 
that in mono- and diethanolamines. The shift in the half-wave potential 
is more complicated due to the following reasons: (1) the degree of irreversi- 
bility increases with the concentration of the amine; (2) there is wave-split- 
ting and (3) the half-wave potentials are almost constant when the solution 
has high concentrations of sodium hydroxide (1-0 M) and a low concentra- 
tion of triethanolamine. But, Lingane’s? equation can only be applied to 
the region where the wave is only slightly irreversible. Under these condi- 
tions, the number of hydroxyl groups calculated for the complex is 4 and 
that of triethanolamine is unity. The formula of the complex can be written 
as [Cd (Trien), (OH),]- ~. 

5. Complex formation in presence of sodium carbonate 

In sodium carbonate solutions also, the cadmium complex with mono- 
and diethanolamine gives a reversible polarographic wave. The formule 
arrived at by applying Lingane’s equation are Cd(Moen),CO, and 
Cd (Dien), COs. 
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When triethanolamine is used as the complexing agent, application of 
Lingane’s equation becomes very difficult, due to the marked interference 
of the second species of the complex and also due to the increasing irreversi- 
bility of the electrode process, with an increase in the conve ntration of tri- 
ethanolamine. In 1-0M carbonate solutions, the interference from the 
second species is not very great. The co-ordination of one molecule of tri- 
ethanolamine with cadmium is obtained by the application of equations for 
irreversible waves. 


6. Complex formation in presence of ammonium chloride-ammonium 
hydroxide 


In presence of ammonium chloride-ammonium hydroxide, the reduction 
process proceeds reversibly with mono- and diethanolamine. The shift 
in the half-wave potential indicates that one monoethanolamine molecule 
has entered into the complex. In the case of diethanolamine the shift is 
too small to suggest the entry of even one molecule of diethanolamine. When 
triethanolamine is used as the complexing agent the polarographic wave 
is reversible only at lower concentrations of the amine. By applying equa- 
tions for irreversible waves? over the concentration range 0-3 to 1:0M of 
the amine we get that one molecule of the amine has entered into the complex, 


7. Thermodynamic data from polarographic measurements 
(i) Hydroxide complex in monoethanolamine—[Cd (Moen), (OH)s]. 
The shift in the half-wave potential of a simple metal ion due to complex 
formation is given by the following expression*:— 
(E s/a)comptex ion —(E1/2)simpte ion = (0°06014/n) log Ka— p (0-06014/n) log C,, 


= E, — E; —p (0-06014/n) log C,, (1) 
where, 


Kg = dissociation constant of the complex, 
C, = excess concentration of the complex forming substance, 
E,. = standard potential of the reaction 
MXp ("P+ +. ne = M(s) + pX-> (2) 
E, = standard potential of the reaction 
M"+ + ne = M(s) (3) 
p and n have the usual significance. 


When the concentration of monoethanolamine and hydroxyl are each equal 
to one molar the equation (1) takes the form 
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(Ex/s)comptex ion (E1/2),smpte fon ~ (0-06014/n) log Kg — E, a E, (4) 
The values of Kg and E, are calculated from the above equation using the 
data given in this paper. (Ey/2)ssmpie ion aS been taken to be — 0-589 V., 
vs. §.C.E. after incorporating theoretical temperature coefficient corrections® 
to the value at 25°C. The standard potential at 30° C. of the reduction 
reaction, Cd**+ + 2e = Cd has been taken to be — 0-402 V.’ (vs. N.H.E). 
The value of 4 F° is calculated from the value of Kg by the following relation 


A F° =—RTInKg (5) 


Cd (Moen), (OH), —> Cd** + 2OH- + 2 Moen 
Kg = 2:244x10-®; 4 F° = 12-00 Kilocals. 


Cd + 2 Moen + 2 OH--—+Cd (Moen), (OH), + 2e; 
Standard potential = + 0-662 V. (vs. N.H.E.). 


(ii) Carbonate complex in monoethanolamine — Cd (Moen), (CQs)). 


Cd (Moen), (CO), — Cd*++ + CO;-- + 2 Moen; 
Kg = 1-:205x10-7; 4 F° = 9-606 Kilocals. 


Cd + 2 Moen + CO, ~—+ Cd (Moen), (COs), + 2e; 
Standard potential = + 0-610 V. (vs. N.H.E.). 
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(iii) Hydroxide complex in diethanolamine — Cd (Dien), (OH),. 






Cd (Dien), (OH), —» Ca** ++ 2 OH- + 2 Dien; 
Kg = 1-225 10-°; 4 F° = 12-38 Kilocals. 
Cd + 2 Dien + 2 OH- —- Cd (Dien), (OH), + 2e; 
Standard potential = + 0-670 V. (vs. N.H.E.). 








(iv) Carbonate complex in diethanolamine — Cd (Dien), (COs), 


Cd( Dien), (CO;); — Cd** + CO, - + 2 Dien 
Kg = 9-55 10-7; 4 F° = 8-359 Kilocals. 


Cd + 2 Dien + CO, -—» Cd (Dien), (CO3); + 2e; 
Standard potential = + 0-583 V. (vs. N.H.E.). 


In the case of triethanolamine it is not possible to calculate the value 
of the thermodynamical constants as the electrode process is mostly irreversi- 
ble. The irreversibility is possibly due to the three hydroxyl groups present 
in triethanolamine which appears to be responsible for the activated nature 
of the reduction process. 
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Fic. 1. Effect of pH 
presence of triethanolamine. 


Curves | to 5 are taken at pH 9-80, 8-65, 7-20, 6-50 and 5- 


The voltage line just before the beginning of the polarogram correspoads to — 0°373 Ve vs. S.C.E. 
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FIG. 2. Effect of sodium hydroxide on the polarographic behaviour of cadmium (1-122 milli- 
molar) in presence of triethanolamine, 


Curves 1, 2 and 3 have 


triethanolamine respectively. 
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the polarographic behaviour of cadmium (1°068 millimolar) in 


been 
Sensitivity: S$/30. 
The voltage line just before the beginning of the polarogram corresponds to -- 0°573 V. vs. S.C.E- 
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taken in 0-1 M sodium hydroxide and 0-30, 0-50 and 0-70 M 
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Fic. 3. Effect of sodium hydroxide (1-0 M) on the polarographic behaviour of cadmium 
(1+122 millimolar) in pre+ence of triethanolamine. 


Curves 1, 2 and 3 are taken with 0-30, 0-50, and 0-70 M triethanolamine respectively. 
Sensitivity : S/30. 


The voltage line just before the beginning of the polarogram corresponds to --0+573 V. 
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SUMMARY 


1. A systematic study of the polarographic behaviour of cadmium has 
been carried out employing mono-, di- and triethanolamines as the com- 
plexing agents at various pH values and in presence of (i) sodium hydroxide, 
(ii) sodium carbonate and (iii) ammonium chloride-ammonium hydroxide. 


2. In base solutions containing triethanolamine, there is a marked 
tendency for the splitting up of the polarographic wave, while this pheno- 
menon is almost absent when mono- and diethanolamines are used as 
complexing agents except when pH is below 11. 


3. The polarographic reduction is reversible in solutions containing 
mono- and diethanolamines and mostly irreversible in triethanolamine solu- 
tions. An explanation has been suggested for the irreversibility of the 
process in triethanolamine solutions. 


4. The formation of the following complexes is indicated from polaro- 
graphic data :— 


[Cd(Moen),(OH)2] [Cd (Moen)(OH).] [Cd (Moen), (COs),] 
[Cd (Moen), (NH;),,]** [Cd (Dien), (OH).] [Cd (Dien), (COs),] 


[Cd (Trien), (OH),]>-- [Cd (Trien),(NHs3)y}** [Cd (Trien), (CO,),.]?-*. 


5. Dissociation constant and the 4 F° value for the dissociation of the 
complex are calculated from polarographic measurements in the case of 
reversible reduction processes. 


6. The use of mono- and diethanolamine in the base solutions is suitable 
for the polarographic estimation of cadmium. 
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ABSTRACT 


The paper contains an analysis of the variation of the midday values 
of foF,, h’F, and hpF, with lunar phase at Ahmedabad during the years 
1954 and 1955 and of foF, alone at Bombay, Madras and Tiruchirapalli 
during 1954. It is found that while the semidiurnal lunar tidal variations 
at Ahmedabad and Bombay agree in phase with those observed at middle 
latitudes, the phase reverses in direction between Bombay and Madras. 
The results are compared with those relating to Huancayo and Singapore. 


INTRODUCTION 


IT is known that lunar tidal variations occur in the upper atmosphere as 
well as in the lower. The amplitudes of variation are much larger aloft and 
the phases are different. Martyn (1949) studied the F2 layer data of a num- 
ber of places in middle latitudes and found that the variation was semidiurnal 
and that the amplitude of variation of foF2 (mean of all hours) range from 
0-06 to 0-142 Mc/s and of hpF2 from 1-7 km to 3-6 km. The maximum 
phases of foF2 and hpF2 were found to be at 10 and 06 lunar hours. 
Appleton and Beynon (1948) using Slough (® = 51° N) data found the ampli- 
tude of lunar variation of foF2 to be 0-05 Mc/s and the maximum phase 
to be at 11 lunar hours, whereas for hpF2, they found corresponding values 
to be 2-0 km and 06 lunar hours. The variation at Huancayo as analysed 
by McNish and Gautier (1949) was different, the amplitude of lunar variation 
of noon foF2 being as high as 0-4 Mc/s (in sunspot minimum years 1941-44) 
with the maximum phase at 04 lunar hours. Thus the lunar variation near 
the geomagnetic equator was apparently higher in amplitude and nearly 
opposite in phase to that at higher latitudes. However, McNish and 
Gautier’s result was obtained by using only noon values of foF2 while 
Martyn used the values at all hours. 04 lunar hours corresponds to nearly 
4 days before the full-moon day. The occurrence of a phase opposition 
near the equator was corroborated by Osborne (1952) from his analysis of 
noon data (1949-51) of Singapore (© = 10-1°S). He found the noon 
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Fic. 1. Variation with lunar phase of midday foF2, hpF2 and h’F2 (10-14 hours) at 
Ahmedabad in 1954 and 1955. 


maxima to occur about four days before full-moon and new-moon days, 
but the amplitude was slightly less than that at Huancayo. Bartels (1950) 
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confirmed the variations at Huancayo. His results were expressed as devia- 
tions from 27-day running averages. He found that the total change in 
foF2 due to lunar tide at noon in southern summer was 1-1 Mc/s while in 
southern winter, it was only 0-2 Mc/s. 


ANALYSIS OF DATA AND RESULTS 


In view of the above results, it was thought worthwhile to analyse the 
F2 data collected at Ahmedabad (® = 13-6°N). The year 1954 was a 
minimum sunspot year with a mean sunspot number 4. As this year was 
almost free from magnetic disturbances, it was considered that 1954 data 
would be appropriate for analysis of lunar daily variation, when other 
influences were at a minimum. To start with, the analysis was made of 
midday foF2, h’f2 and hpF2, taken to be the mean of the values at 10-14 
hours. The days in all the months having the same phase (x) of the moon 
were grouped together and 24 such groups were obtained for the different 
phases of the moon. For each phase, the mean midday values (mean of the 
5-hourly values at 10-14 hours) for the whole year was found. When- 
ever there was loss of record or blanketing by Es or failure of equipment, 
reasonable interpolation was made for the missing value. All the days for 
which data could be had were taken into account and seasonal changes were 
eliminated by taking yearly means. The results are given in Fig. 1, which 
clearly shows that the maximum amplitude occurs about 2 days after the 
full-moon and new-moon days for foF2 and about 5 days after full or new 
moon days for h’F2 and hpF2. The results are expressed as percentage 
deviations of 3-day mean centred at each lunar day from the 15-day running 
averages. Some points lie outside the curve and these are found to be just 
those for which the number of observations were few. In terms of lunar 
time, the phase of foF2 comes out to be 10 lunar hours and that for h’F2 
and hpF2 about 07 lunar hours. The corresponding amplitudes are found 
to be about 4-5% in foF2 and about 1-7% in heights. The results for the 
year 1955 on the rising part of the solar cycle are shown in the same diagram. 
The maximum in 1955 has shifted earlier by a day compared to 1954, both 
in foF2 and h’F2. Fig. 2 gives the mean curves of lunar variation in foF2 
and h’F2 for both the years together. The points for each year are marked 
differently. It will be seen that the phase agrees well with those observed 
at higher latitudes by Martyn, but the amplitude is larger. One reason for 
this may be that we have considered midday values only. 


It was thought that it would be interesting to examine the lunar variation 
in foF2 at other Indian stations nearer the geomagnetic equator, viz., Bombay, 
Madras and Tiruchirapalli. The midday values of foF2 for the year 1954 
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Fic. 2. Variation with lunar phase of midday foF2 and h’F2 (10-14 hours) at Ahmedabad 
in 1954-1955. 

foF2: ---- 1954, #444 1955 ; h’F2: ooco 1954, xxxx 1955. 

Full line curve is the mean for 1954-55. 
at the above stations collected by All-India Radio were analysed in the same 
manner, and the results are given in Fig. 3. It is interesting to note that 
the phase at Bombay (® = 9-5° N) is practically the same as at Ahmedabad, 
whereas at Madras (® = 3-1° N), the phases are almost opposite. This 
shows that north of equator, a phase-reversal takes place somewhere between 
9-5°N and 3-1° N geomagnetic latitude. Tiruchirapalli which is still near 
the equator has a maximum about 24 days before full or new moon as against 
about 34 days at Madras. Another point is that instead of an increase in 
amplitude towards the equator as observed by others, we find that the 
amplitude decreases after reaching a maximum. Observations at Ahmedabad 
were made with an automatic ionospheric recorder whereas at the other 
places they were made manually. The accuracy of reading in foF2 was 
0-1 Mc/s at all the places. 


It is interesting that although Bombay has nearly the same geomagnetic 
latitude to the north of the geomagnetic equator as Singapore is to the 
south, there is a difference between the lunar variations at the two places. 
While the equatorial type of variation is observed at Singapore, the middle 
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Fic. 3. Variation with lunar phase of midday foF2 (10-14 hours) at different latitudes in 
1954. A: Ahmedabad, B: Bombay, M: Madras, T: Tiruchirapalli. 




















latitude type occurs at Bombay. The change to the equatorial type occurs 
between Bombay and Madras, perhaps nearer Bombay than Madras. 
Singapore is however only 1° 19’ north of the geographic equator. 


To facilitate ready reference, the results of the analysis of the second 
harmonic components for 1954 are collected in Table I below; those of 
Singapore relate to 1948-51 and of Huancayo to 1941-44. 
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TABLE I 
Amplitude P2 and Phase t2 of lunar Variation in foF2 





foF2 hpF2 





Place Gceomeg. 
lat. 100P2 P2 t? lunar P2 t2 lunar 
mean foF2 Mc/s hours km hours 





Ahmedabad .. 13-6°N 3-9 0-35 10-2 3-0 7-2 


Bombay io IN 2-6 0-24 10-0 
Madras . cr 1-8 0-14 3°8 
Tiruchy ui or 1-8 0-13 2°8 
Huancayo .. 0-6°S 5-0 0-40 4-0 
Singapore .. I9°s (3-0) (0-30) (4-0) 





The non-agreement of the results of Bombay and Singapore suggests 
that the lunar variation does not depend on geomagnetic latitude only. If 
we consider two stations according to geographical latitudes, Bombay 
(geogr. lat. 19° N) falls in the middle latitude group while Singapore 
(geogr. lat. 1-3° N) is in the equatorial group. 

Bartels has shown that at Huancayo, the total change of noon foF2 
during a lunar cycle in southern summer is nearly 5 times that in southern 
winter and Burkard (1951) has shown that the semidiurnal lunar variation 
of foF2 is pronounced during the daytime only and practically absent during 
the night. The significance of the phase-reversal of the lunar semidiurnal 
variation between the latitudes of Madras and Bombay and at some Iati- 
tude to the south of Singapore remains to be explained. An analysis 
of the available Indian data to study the lunar variations in different seasons 
of the year is being undertaken. 
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